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The phase transition to superfluidity and the BCS-BEC crossover for an ultracold gas of fermionic 
atoms is discussed within a functional renormalization group approach. Non-perturbative flow 
equations, based on an exact renormalization group equation, describe the scale dependence of 
the flowing or average action. They interpolate continuously from the microphysics at atomic or 
molecular distance scales to the macroscopic physics at much larger length scales, as given by the 
interparticle distance, the correlation length, or the size of the experimental probe. We discuss the 
phase diagram as a function of the scattering length and the temperature and compute the gap, 
the correlation length and the scattering length for molecules. Close to the critical temperature, we 
find the expected universal behavior. Our approach allows for a description of the few-body physics 
(scattering and molecular binding) and the many-body physics within the same formalism. 



PACS numbers: 03.75.Ss; 05.30.Fk 



I. INTRODUCTION 

The binding energy of molecules formed from two 
fermionic atoms may depend on external parameters such 
as a magnetic field B. We discuss a gas of two degen- 
erate species of atoms (typically hypcrfine states), where 
bound molecules form for B < Bq, whereas for B > Bq 
only unstable resonance states are present. At the Fes- 
hbach resonance B = Bo, the scattering length a di- 
verges. At low temperatures T and for large B — B , 
the molecules play no role and one expects superfluid- 
ity of the BCS type [J. In the opposite case of large 
negative B — B , the bosonic molecules dominate at low 
temperatures T and Bosc-Einstcin condensation (BEC) 
should occur. In the vicinity of the Fcshbach reso- 
nance, a continuous crossover from BCS to BEC con- 
densation has been predicted Q. Several recent experi- 
ments have established this overall picture In paral- 
lel, a considerable theoretical effort has been devoted to 
a quantitatively precise theoretical understanding of this 
phenomenon [5l424|. 

A major motivation for the theoretical effort is due to 
the possibility of using this system as a testing ground for 
non-perturbative methods in many-body theory, which 
go beyond the mean-field approximation. Quantitatively 
reliable methods for complex many-body systems are cru- 
cial for many areas in physics, from condensed matter 
phenomena to nuclear or elementary particle physics. 
For a non-rclativistic system, as in our case, the results 
of such methods may be tested, in principle, by a huge 
number of experiments and observations. However, for 
a quantitatively precise test the underlying microphysics 
must be known with sufficient precision. This is often 
not the case for real physical systems in condensed mat- 
ter physics. A notable exception is the universal physics 
near a second order phase transition, where the criti- 
cal exponents and amplitude ratios no longer depend on 



microscopic details. The quantitative computations of 
these universal quantities by the renormalization group 
has been one of the major breakthroughs in theoretical 
physics [H,[26[. 

In cold atom gases, the microscopic physics is well un- 
der control, since it can be directly related to atomic and 
molecular properties and measured by scattering experi- 
ments. This gives rise to the hope that experimental tests 
become now available for the whole connection between 
microphysics on the one and the macroscopic quantities 
as thermodynamics, condensation phenomena, the cor- 
relation length, defects etc. on the other side. Further- 
more, the effective strength of interaction can be tuned 
experimentally, which offers the fascinating perspective 
of following the tests of theoretical methods from the 
well-controlled perturbative domain for small scattering 
length to the difficult non-pcrturbativc region for large 
scattering length. 

First results for the BCS-BEC crossover are now avail- 
able for a number of theoretical methods that go be- 
yond mean-field theory. Quantitative understanding of 
the crossover at and near the resonance has been devel- 
oped through numerical calculations using various quan- 
tum Monte-Carlo (QMC) methods 0-11]. Computations 
of the complete phase diagram have been performed from 
functional field-theoretical techniques, in particular from 
e-expansion 9l-[l2^. 1/iV-expansion [l3|, [l4|, t- matrix ap- 



proaches [15] . Dyson-Schwin ger equations [lg, 117] . 2- 
Particle Irreducible methods jl8j |. and renormalization- 
group flow equations fl"9U24l ] . 

In this paper, we give an account of the func- 
tional renormalization group approach for the BCS-BEC 
crossover which is based on the flowing action (or aver- 
age action) (27]]. We concentrate here on the detailed 
description of the formalism which underlies the results 
presented in ref. [2(1 [23|- Furthermore, we extend 
the results of these papers by results for the correlation 
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length and the self-interaction of the bosonic molecules, 
including the behavior near the phase transition. 

Our paper is organized as follows. Sect. Ql] focuses 
on the method which we use in order to tackle the 
crossover problem, the functional renormalization group. 
Our starting point is an exact functional renormalization 
group equation for the average potential [27j . We specify 
the approximation scheme used for our practical compu- 
tations and present the non-pcrturbativc flow equations. 
Details and various quantities needed for a numerical so- 
lution can be found in a series of appendices A-F. This 
section provides the technical basis for our work. 

The physics of the crossover problem is addressed in 
the following sections. It can be structured into two 
parts: Few-body and many-body problem. In our for- 
malism, both parts find a unified description. In Sect. 
IIIII we discuss the phase transition to superfluidity at 
low T within our formalism. We present results for the 
temperature dependence of the gap and the correlation 
length for T < T c . For a comparison with experiment, 
we need to relate the microscopic parameters appearing 
in the microscopic action to physical observables as the 
scattering length or the molecular binding energy. The 
latter quantities concern the behavior of a two-atom sys- 
tem. In our formalism, this can be described as the vac- 
uum limit where density and temperature approach zero. 
This is discussed in Sect. II VI Sect.lVlextends the simplest 
truncation (discussed in Sect. [TTJ) in order to include the 
particle- hole fluctuations (23|. In Sect. I VII we demon- 
strate the capacity of the model to incorporate directly 
the universal critical physics at the phase transition. We 
discuss the behavior of the scattering length for molecules 
as the critical temperature is approached. In Sect. IVIII 
we address the question of a quantitatively precise com- 
putation of the density. The density sets the only scale in 
the unitarity limit where the scattering length for atoms 
diverges. It enters directly the theoretical predictions 
as the ratio T c /Tp, where the Fermi temperature T-p is 
given by the density. The determination of the density by 
a direct computation of the thermodynamic potential on 
the chemical potential fi leads to a lower value than the 
"free particle density". As a result, the value of T c /Tp is 
somewhat higher. We present conclusions in Sect. I Villi 

II. FUNCTIONAL RENORMALIZATION 

A. Effective Action and Potential 

We express our microscopic model in terms of the ac- 
tion 

S[tp, 1>] = ^ dr J d 3 x(V (<9 T - A - n) V> 
+<p* (d T - + v — 2/zj tp 
-^(^e^-^er)), (1) 



where ip T ev(j = 2-0i?/>2, ip^eip* = 2^*^>| and A is the 
Laplacian. This model describes nonrelativistic fermions 
in two different states, a composite bosonic field and a 
Yukawa type interaction between them. We represent 
the fermion field by the complex two-component spinor 
ip = (tpi,tp2) (e is the totally antisymmetric symbol) 
and the bosonic complex scalar field by tp. Depending 
on the situation, if can describe a bound state of two 
fermions (molecule), the closed channel state of a Fesh- 
bach resonance or a Cooper-pair like collective state of 
two fermions. As long as its dynamics plays no role, it 
can also be considered as an auxiliary field that is used 
as an effective parameterization of a purely fermionic in- 
teraction. 

We work with the formalism of quantum statistics in 
the grand canonical ensemble with chemical potential /z. 
In addition to the spatial position x, the fields depend 
on Euclidean time r. At nonzero temperature, this Eu- 
clidean time is wrapped around a torus of circumference 
/? = j!. This implies for the boson field <p(x, (3) = <p(x, 0), 
while the fermion field gets an additional minus sign due 
to its Grassmann property ip(x, (3) = —i/)(x, 0). Here, and 
in the following we choose our units such that h = ks = 1- 
In addition, we also rescale time, and the fields in units of 
the fermion mass M such that effectively 2M = 1. App. 
A gives a brief summary how eV units for the physical 
observables are recovered from our units. 

At tree level, the Yukawa term leads to an effective 
four fermion interaction ipi4>2 yi y2i with coupling 
strength 

A^efi = — -T — nr^r ( 2 ) 

v - 2[i + iq + -|<r 

Here, qo and q denote the center of mass Euclidean fre- 
quency and momentum of the scattering particles ipi and 
ip%. The real time frequency u> obtains by analytic contin- 
uation, iqo — > — uj, such that the microscopic dispersion 
relation for the molecules is ui = q^ /2 + v — 2/i. In the 
limit fit — > oo,^ — > oo,hz,/v — > const, one can neglect 
go j Q 2 and \i in Eq. ([1]), such that the effective fermionic 
interaction becomes pointlikc, cf. [l6| and Sect. IIV CI 

Our aim is the computation of the quantum effec- 
tive action T[<p,?p]. It has a structure similar to the 
microscopic action S, but with microscopic couplings 
and fields replaced by renormalized macroscopic cou- 
plings and fields. Furthermore, the couplings will depend 
on momentum, and T contains additional couplings not 
present in S. The quantum effective action generates 
the one-particle-irreducible vertices (amputated scatter- 
ing vertices) and (inverse) full propagators. The field 
equations derived from an extremum principle for T are 
exact equations - all contributions from quantum and 
thermal fluctuations are included. All scattering infor- 
mation can be easily extracted from T as well as all ther- 
modynamic relations. 

For constant scalar fields <p(x, r) = <p and tp(x, t) = 0, 
the effective action reduces to the effective potential 
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U(ip), since T = (V 3 /T)U(ip) with V 3 being the vol- 
ume. The U(l) symmetry of global phase rotations im- 
plies that U can only depend on the invariant p = ip*(p. 
The field equation dU /dp = determines the location 
po of the minimum of U. Any nonzero tp amounts to 
a spontaneous breaking of the £7(1) symmetry. The as- 
sociated Goldstone boson is responsible for superfluidity. 
Therefore, ipo is an order parameter - more precisely 
2po = 2pqPq is the supcrfluid density. The system is 
in the superfluid phase if p (T,p) > 0, and in the disor- 
dered or symmetric phase if po{T, p) = 0. The knowledge 
of U(p; T, p) thus permits a straightforward extraction of 
the phase diagram. Furthermore, the value of U at the 
minimum corresponds to the pressure 




v{T,p) = -U(p \T,p). 



(3) 



(We choose a free additive constant in U such that 
p(0, 0) = 0.) The grand potential or Landau thermody- 
namic potential $g is related to U(po) by 



$G = V 3 U{po), 



(4) 



and we obtain the usual relations for the particle density 
n, entropy density s and energy density e: 



dp 
dp 
dp 
df 



dU 
dU 



-p + Ts + pn = U (pa) + Ts + pn. 



B. Flow equation for the average potential 



(5) 



FIG. 1: Flow of the effective action in the space of couplings. 
The flow starts at an initial scale A, where all fluctuations 
are suppressed and the effective action reduces to the classi- 
cal action. By solving the FRG equations we follow the flow 
in parameter space down to k = 0. Two-body scattering ob- 
servables are measured in the physical vacuum defined as the 
state with n = T = 0, but at k = - all vacuum fluctuations 
are included. The classical action at k — A is chosen such 
as to match the two-body observables in the physical vacuum 
state (dashed line). We then switch on temperature and den- 
sity to solve the many-body problem. The many-body flow 
deviates from the vacuum flow when the cutoff wavelength 
becomes comparable to the characteristic thermodynamic 
wavelengths, i.e. the de Broglie wavelength A c ib ~ T~ 1//2 and 



the mean interparticle spacing d > 



-1/3 



(solid line). 



The dependence of Uk on the "average scale" k obeys 
an exact renormalization group equation or flow equation 
[H, with q= (uj,q) 



The effective potential U contains all the essential in- 
formation for the description of the equilibrium proper- 
ties of a homogeneous system. Its computation is, how- 
ever, a challenge, since quantum and thermal fluctuations 
on all momentum scales have to be computed. The av- 
erage potential Uk includes only the effects of bosonic 
fluctuations with momenta q 1 > k 2 and fermionic fluctu- 
ations with {cp 2 — Py \ > k 2 , withpp the Fermi momentum. 
It is therefore a type of coarse-grained free energy. If k is 
large enough, the complicated infrared physics for bosons 
and physics close to the Fermi surface for fermions is not 
yet included, and one expects that the complexity of the 
problem is reduced. In the limit k — > oo, no fluctua- 
tion effects are included at all, and the average potential 
becomes the classical potential extracted from (fT]), 

U k ^oo{p) = U c i(p) = (v - 2p)p, p = ip*ip. (6) 

On the other hand, for k — > all fluctuation effects are 
included and Uk=o equals the effective potential U. The 
challenge therefore is to construct an interpolation from 
the microscopic potential Uk^oo to the effective potential 



-tT il [G^( l p;q) + R^\q)]- 1 d k R^\ q )}. 



(7) 



The first piece accounts for the bosonic fluctuation and 
the second for the fermions. Eq. ([7]) is an exact functional 
differential equation. It includes all orders in a pcrturba- 
tive expansion as well as all non-perturbative effects. 

In order to understand the structure of this equa- 
tion we first include only the bosonic contribution. Our 
model reduces in this case to a gas of bosons, for which 
the flow equations [28| and thermodynamics [29| have 
been discussed extensively in presence of a repulsive in- 
teraction and for arbitrary dimension. In Eq. ([7]) we 
encounter the propagator G v of the scalar field in an 
arbitrary background ip. We use a two-component ba- 
sis of real fields ipi,if2, related to the complex field by 
tp(x) = + i(p 2 (x)). Correspondingly, G~ x is a 

2x2 matrix and tr v denotes the trace in this space. In 
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the classical approximation, the inverse propagator reads 



q 2 /2 + V-2/J, 



-So 



<?0 



q 2 /2 + v-2fi 



(8) 



In the absence of an interaction between the bosons, it 
is independent of (p. For v = 2/U, the classical propa- 
gator G VtC i becomes singular for q = (go, q) — > 0. It is 
regularized by adding a momentum-dependent infrared 
cutoff R k (q) (a unit matrix in ipi,ip2 space) which has 



the general properties 
lim R k (q) ' 



lim R k (q)^0. 

i 

(<eh 



(9) 



For k > the combination (G" 1 T 1 

remains finite 

even for q z /2 + 1/ — 2/i = 0. In principle, Eq. (0 holds for 
arbitrary cutoff functions Rk(q) that fulfill Eq. ©. For 
approximate solutions of Eq. an appropriate choice is 
important, however, and we discuss our choice of Rk in 
App. O 

The structure of the flow equation already becomes 
clear in a perturbative one loop truncation: if we replace 
G v by G v ,d in the flow equation, the inverse propagator 
is independent of k such that the k integration can be 
performed analytically 



U k = U d + 



tr^HG" 1 +R k ] 



(10) 



With j. = / A7(20 3 , J qo = T£„, 90 = 27rnT, we 
recognize the regularized one-loop contribution. For k = 
and v = 2fi, one rediscovers the pressure of a gas of 
free bosons. For T < T c , this describes Bose-Einstein 
condensation. 

The exact flow equation is obtained if we perform a 
"rcnormalization group improvement" of the propagator 
and replace in Eq. (0 the classical inverse propagator by 
the full inverse propagator in the presence of a "back- 
ground field" ip. The full inverse propagator is given by 
the second functional derivative of the average action Tk 
with respect to the background fields <p(x). This average 
action generalizes the effective potential in the presence 
of derivative terms for the fields. We will discuss this in 
more detail below. What is important at this stage is 
that G" 1 typically depends on ip, thus generating a non- 
trivial ^9-dependence in the flow equation ([?])• Indeed, in 
the presence of effective bosonic interactions Tk will not 
be quadratic in (p anymore. 

If one neglected the fc-dependence of G" 1 , Eq. (fTU|) 
would give the one-loop perturbative result for interact- 
ing bosons. However, beyond this approximation G^ 1 
will depend on k since only fluctuations with momenta 
<f > k 2 are included in its computation. Then Eq. (J7J) 
becomes a non-linear differential equation that has to be 
solved numerically. Nevertheless, its structure always re- 
mains a one-loop equation with only one momentum inte- 
gration. The only difference from one loop perturbation 



theory arises from the use of the full field and momen- 
tum dependent propagator. Higher (perturbative) loops 
as well as non-perturbative effects are created due to the 
non-linear nature of the flow equation. 

In the BEC-limit the fermion chemical potential \i as- 
sumes negative values, as appropriate for a gas of bosons. 
Starting from the action in Eq. ([1} we could perform the 
Gaussian functional integral for the fermions, resulting 
in an effective action involving only the bosonic field ip. 
In general, the effective bosonic action is a complicated 
non-local object. However, for large enough negative /x 
(as compared to the typical frequencies u and kinetic 
energies q 1 /2) the bosonic action becomes effectively lo- 
cal. We could then use the purely bosonic flow equa- 
tion, starting with a non-polynomial classical potential 
U\(ip). We therefore expect to recover in the BEC-limit 
many features of a boson gas with repulsive interactions. 
(There may be quantitative differences as compared to 
the case of a pointlike interaction between two bosons, 
due to the non-polynomial U (ip) which also depends on 

Away from the BEC-limit the fermion fluctuations will 
play a crucial role - they dominate in the BCS-limit. We 
therefore have to solve the full Eq. ([7]). The fermionic 
contribution in Eq. (J7J is of a similar nature as the 
bosonic contribution, with the characteristic minus sign 
for a fermion loop and G^ 1 being the full fermionic in- 
verse propagator. The fermionic cutoff R k is chosen 
such that it regularizes the momenta close to the Fermi 
surface rather than the small momenta, see App. [C] 

The BCS limit obtains for fx > and large positive v 
and hfp/v. The qualitative features can now be under- 
stood in the limit where the boson fluctuations can be 
neglected. The inverse fermion propagator G^ 1 will con- 
tain a piece ~ h v ip due to the Yukawa interaction in Eq. 
(p}. This induces contributions to the effective potential 
U(tp), which lead to a minimum away from tp = if the 
temperature is low enough. In the classical approxima- 
tion the fermion inverse propagator reads 



-h 



v eip 

q 2 - a* 



m - (q 2 - m) 

h v tip 



(11) 



This is a 4 x 4 matrix, where we have suppressed the 
spinor indices. (The symbol e a p acts on the spinor indices 
a = 1,2.) In the limit where the /c-dependence of h v is 
neglected we can again integrate in order to obtain the 
one loop result for k = 



1 

{v A - 2n)ip*ip 



tr,/, In G 



(12) 



2T / In 



cosh 



(yV - A*) 2 + K^P/2T 



Condensation occurs if the minimum of Uq occurs for 
\ifo\ > 0. This always happens for small enough temper- 
ature T . Determining from Eq. (|11[) the critical temper- 
ature T c for the onset of a condensate ipo ^ yields the 
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BCS value. (For this purpose we relate va to the scat- 
tering length a in Sect. [TV]) For quantitative accuracy 
the bosonic fluctuations cannot be neglected even deep in 
the BCS region. We will see in Sect. [V] that they induce 
in our formalism the corrections from the particle-hole 
fluctuations. 

From Eq. (fT2")) we can read off the fermion dispersion 
Eq = J(JP- M) 2 + ffypfrfo- Thus, A = hptpo repre- 
sents the U(l) symmetry breaking mass term or gap for 
the single particle excitations. In our framework, it is 
proportional to the bosonic field expectation value, such 
that spontaneous symmetry breaking po — Vo'Po 7^ 
and a gap for the fermions coincide. The actual values 
of h v , po are, however, determined beyond the mean field 
approximation, cf. e.g. Eq. (j2"Tj) . 



C. Truncation 

The full propagators G v ^ remain complicated objects. 
It is at this level that we need to proceed with an ap- 
proximation by making an ansatz for their general form. 
This amounts to a truncation of the most general form 
of the average action r&, which is the generalization of 
the average potential Uk to a functional of space- and 
time-dependent fields tp(x),ip(x). In other words, is 
the quantity corresponding to the effective action T in 
the presence of the infrared cutoff Rk, which suppresses 
the fluctuation effects arising from momenta q 2 < k 2 . 

A truncation of Tk should respect the symmetries, 
which we discuss in App. [Bj Since fluctuations for the 
derivative terms will turn out to be important for the 
composite bosonic field we have to distinguish carefully 
between microscopic (bare) and renormalizcd fields and 
couplings. We use in this paper the lowest order in a 
derivative expansion which, in terms of bare quantities, 
reads 



r fc = 



i 



^\u k (p) + {p*(Z v d T -^A v &)y 



(13) 



The bar denotes bare fields and couplings, which can 
be related to renormalized quantities by absorbing the 
running gradient coefficient A v by a field rcscaling ip = 
l li 

A v <p (wavefunction rcnormalization) . As we neglect 
the fermionic wave function rcnormalization in this work, 
there is no need to distinguish between bare and renor- 
malized fermion fields here. In terms of renormalized 
fields, the truncation then reads 



Tk = 



{u k { P ) + v *{s v d T 

X 

+V f (<9 T - A — n)tp — ^(^ilFeil) - tprfl 



with the renormalized couplings S v 
htp/ yAZ, order parameter p = A v p, and effective po- 



^tp/Ap, 



h v = 



tcntial U(p) = U(p). In addition to the average po- 
tential Uk, the wave function renormalizations Z V ,A V 
and the Yukawa coupling h v depend on k. In this ap- 
proximation the renormalized inverse boson propagator 



G„ 



G /A v reads 



G~ 



W k ' + q 



2 /2 



Here, primes denote derivatives with respect to p. The 
fermionic inverse propagator is given in Eq. (Jill) , where 
we use a basis (ip, ip*). 

We emphasize that Uk appears on the right-hand side 
of Eq. ([7]): at any given scale k the flow is determined 
by the average potential Uk at the same scale k rather 
than by the classical potential. The same holds for the 
/c-dependent Yukawa coupling h v . The right-hand side 
of the flow equation is only sensitive to the physics at the 
scale k. In particular, the structure of the flow equations 
implies that the momentum integrals in Eq. (UJ are dom- 
inated by a small range q 2 rj k 2 . This is the main reason 
why a derivative expansion often works rather well. In 
principle, A v , Z v and h v are momentum-dependent func- 
tions. If this momentum dependence is not too strong, 
and the effective range of momenta in the integral ([7]) 
remains rather small, a momentum-independent approx- 
imation is a good starting point. The truncation Eq. (| 14[) 
forms the basis for the results presented in this paper. 
A similar truncation has been advocated for attractive 
fermions on the lattice in Ref. [3Cj . In Sect. |V]we extend 
it to assess the effect of particle-hole fluctuations. 

The truncation in Eq. (|13[) has been extended in ref. 
[3l| where renormalization effects on the fermionic self- 
energy and in particular a wavefunction rcnormalization 
factor for the fermions have been included. For a detailed 
discussion we refer to ref. 1311] . 



D. Non-perturbative flow equation 

For our choice of Rk and with the approximation (1151) . 
we can perform the momentum integration and the Mat- 
subara sums explicitly 



kd k U k = VA v pU' k + 



V2k 5 
3tt 2 S u 



(1-2^/5) 



.(0) 



-—l(u)s (0) 



(16) 



m = (0(m + i)(a + i) 3/2 -^-i)(A-i) 3/2 ). 

Here we use the dimensionless chemical potential p = 
(14) p/k 2 , the ratio S v = Z v jA^ and the anomalous dimen- 



sion 77^4^ = — dlnA v /dhik. The threshold functions sb 
and sf depend on Wi = U' k /k 2 , w 2 = (U' k + 2pU k )/k 2 , 
W3 = h^p/k 4 , as well as on the dimensionless tempera- 
ture T = T/k 2 . They describe the decoupling of modes 
if the effective "masses" or "gaps" Wj get large. They are 
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o 



»- + dt 



•»--•»'--»- 



d t >-»-- = 



FIG. 2: Diagrammatic representation of the flow equations 
for bare quantities in the symmetric phase. The flow of the 
inverse propagator is displayed in the first line. The frequency 
and momentum dependence of the loops is projected accord- 
ing to Eqs. (|24|) . The fermionic graph, present in the mean- 
field treatment taking Gaussian bosonic fluctuations into ac- 
count, is supplemented with a feedback of bosonic fluctua- 
tions due to the second graph. These effects are important to 
capture the expected long-wavelength behavior close to the 
second-order phase transition, whereas this graph vanishes 
in the physical vacuum n — T — 0. The Yukawa vertex in 
the second line does not flow within our truncation due the 
U(l) symmetry. This changes in extended truncations [23| . 
Moreover, in the symmetry-broken phase the flow equations 
become more complex, due to processes describing scattering 
off the condensate. In this case, also the Yukawa vertex ac- 
quires a nonzero flow. The dt derivative acts on the cutoff 
functions in the regularized propagators only. 



normalized to unity for vanishing arguments and T — > 
and read 



,(0) 



,(()) 



1 + W\ 



w 2 



1 + U> 2 
1 + W! 



+ N B (^/l + w 1 x/l + w 2 /S ip ) 
I - N F {Vl + w 3 ) 



(17) 



(For Sg , only all its p derivatives vanish for w\ ~ w 2 — > 
oo. The remaining constant part is a shortcoming of the 
particular choice of the cutoff acting only on spacelike 
momenta.) 

In Eq. (fT7| the temperature dependence arises through 
the Bose and Fermi functions 



N B/F (e) 



1 



e/T T X 



(18) 



For T — > the '"thermal parts" ~ A^f vanish, whereas 
for large T one has 



„(°) 



2fS v {l+wi)- l {l + w 2 y 1 . (19) 



In this high-temperature limit the fermionic fluctuations 
become unimportant. For the boson fluctuations only the 



n = Matsubara frequency contributes substantially. In- 
serting Eq. (|19p into Eq. (JTS]) yields the well known flow 
equations for the classical three-dimensional scalar the- 
ory with U(l) symmetry [27], [H[. In App. |D]we derive 
the flow equations (fT6f and discuss the threshold func- 
tions sb and sf more explicitly. 

We will further truncate the general form of the aver- 
age potential 



-p k + m v (p-po 
-n k 5p + a k {p - p )6p 



Po) 2 



(20) 



Here m 2 „ 



\{k), \ v = \ v {k), p = po(k) as well as 



Pk, ctk depend on k. All couplings depend also on T 
and a reference chemical potential p via p = p a + 5 p. 
(The differentiation with respect to p does not act on 
Po but rather on Sp, which we set to zero after the p 
differentiations.) In the symmetric regime, we have po = 
0, whereas the location of the minimum po(k) becomes 
one of the parameters in the superfluid regime (po(k) > 
0). In the latter case, we have to take = such that 
Po(k) indeed corresponds to the minimum of Uk for given 
p. We recall that Pk^o is th e pressure. 

An expansion of the effective potential around the min- 
imum as in Eq. (|20[) is expected to work reasonable for 
the description of second order phase transitions. In more 
general situations where also first order phase transitions 
play a role, one might use a truncation where the full 
function Uk(p) is kept. The partial differential equation 
([7]) can then be solved numerically, for example using a 
discretization technique. In the present work we concen- 
trate on the description of the second order phase transi- 
tion of spin-balanced Fermi gases so that the expansion 
in Eq. (|2T)|) is sufficient. 

In the truncation of the effective potential in Eq. (|20l) , 
one might include as a next step a term ~ (p — p ) 3 - 
We have not checked the influence of such a term for the 
BCS-BEC crossover. For a Bose gas investigated with 
a similar truncation [28| , the influence of such a term is 
quite modest, however. 

The flow equations for pk,m/^ or po(k), and are 
given by 



dkPk 



-d k U k 



d k ml = d k U' k 



PQ' 
k I p=0 



d k p 
d k \ v 



for po = 0, 
for po > 0, 



(21) 



Taking a derivative of Eq. (fl6|) with respect to p one 
obtains for T = 



kd k U' k = r ]A jU k -pU' k r ) + 



V2k 



1 



d + 2 



o) 



MUZ) 2 (4 1 , 



3s 



(0,1) 
B.Q 



"7t(3)J0,1) 



B.Q 



F,Q- 



(22) 
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The threshold functions s^q , Sgo i an< ^ s fq are defined 
in App. [D]and describe again the decoupling of the heavy 

modes. They can be obtained from p derivatives of 



and 



JO) 



Setting p = and T — > 0, we can immediately 



infer from Eq. 
regime. 



the running of nr^ in the symmetric 



v - kd k U' k = T] Aip m v + 7^h v l{p) 4,q(^3 = 0). 



kdkm 

(23) 

One can see from Eq. (f2"B"| that fermionic fluctuations 
lead to a strong renormalization of the bosonic "mass 
term" m 2 . In the course of the renormalization group 
flow from large scale parameters k (ultraviolet) to small 
k (infrared) the parameter m 2 ^ decreases strongly. When 
it becomes zero at some scale k > the flow enters the 
regime where the minimum of the effective potential Uk 
is at some nonzero value po- This is directly related to 
spontaneous breaking of the U(l) symmetry and to lo- 
cal order. If po ^ persists for k — > this indicates 
superfluidity. 

For given A v , Z v ,h v , Eq. (|16[) is a nonlinear differ- 
ential equation for Uf~, which depends on two variables 
k and p. It has to be supplemented by flow equations 
for Ap, Z v , h v . The flow equations for the wave func- 
tion renormalization Z v and the gradient coefficient A v 
cannot be extracted from the effective potential, but are 
obtained from the following projection prescriptions, 

d 



d t z v = -9 t _(P v ) 12 ( go ,o)| go=0 , 



d t A 



2d t ^r(P v ) 22 (0,q)L 



(f=0' 



(24) 



where the momentum dependent part of the propagator, 
P v {q) = G-\q) - G-\0) is defined by 



8tp a (q)6tp b (q') 



(P v )ab(q)S(q- 



(25) 



The computation of the flow of the gradient coefficient is 
rather involved, since the loop depends on terms of differ- 
ent type, ~ (q ■ p) 2 , q 2 , where p is the loop momentum. 
The computation is outlined in App. El together with 
the calculation of Z v and the respective flow equations. 
In App. [F]we describe the flow equation for the Yukawa 
coupling h v . 

We provide the diagrammatic representation of the 
flow equations considered here in the symmetric phase 
in Figs. ECU 



E. Scattering length and concentration 

Besides the thermodynamic parameters temperature T 
and the atom density n, our model contains free param- 
eters which characterize the microscopic action T\ and 
constitute the initial values for the flow. The most im- 
portant one is the "detuning" z/, which can be related to 



\ & 






d t X = 


< 


\ + &t ^ x 




> 





FIG. 3: Diagrammatic representation of the flow equation 
for the bosonic self-interaction in the symmetric phase. As 
in Fig. [5] we observe a feedback due to bosonic fluctua- 
tions. They have an important impact on the ratio of bosonic 
to fermionic scattering length as described in Sect. IV G 
extracted in the physical vacuum, as well as on the long- 
wavelength physics close to the phase transition. 



the scattering length a for atom-atom scattering. The 
relation between v and a only concerns the two-body 
problem. We will redcrive the well known function a(v) 
within our formalism in Sect. IIV1 Relations between 
macroscopic quantities only involve dimensionlcss quan- 
tities, once appropriate units are set. We will therefore 
use instead of the scattering length the dimensionless 
concentration c = etfcp, where we define the Fermi mo- 
mentum kp in terms of the density, n = fcp/(37r 2 ). (This 
is a formal definition that we also use for T ^= and in 
the BEC regime where no Fermi surface exists.) For a 
small value of \c\ the scattering length is small as com- 
pared to the inter-particle distance and the system is 
weakly interacting. In contrast for \c\ > 1 we have to 
deal with a strongly interacting system. The "unitarity 
limit" corresponds to diverging a or = 0. Positive 
c corresponds to the BEC side of the crossover where 
stable molecules exist in vacuum. The BCS side of the 
crossover, where the condensation phenomena are related 
to "Cooper pairs" of atoms, is described by negative c. In 
summary, two of the parameters appearing in T\, namely 
v and /i, will be replaced by a and &f, and chosen in order 
to reflect a given concentration c. 

The microscopic model contains the Yukawa or Fesh- 
bach coupling as an additional free parameter. We will 
consider the limit of a strong Yukawa coupling, which 
corresponds to a broad Fcshbach resonance. In this limit 
the flowing Yukawa coupling approaches rapidly a par- 
tial fixed point [2lj , such that its precise microsopic value 
is not relevant. For a broad Feshbach resonance, many 
other possible microscopic parameters play no role either. 
This aspect of universality of broad Feshbach resonances, 
which also holds away from the unitrarity limit, has been 
extensively studied in [2l| . 

To end this section let us make some remarks on the 
ultraviolet scale A. For our numerical solution of the flow 
equation we have to choose some value of A but as long 
as it is large enough, its precise value is not important. 
The reason is that for infrared cutoff scales k larger than 
the scales given by the temperature k 2 > T and the den- 
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sity k 3 > n, the flow equations have the same form as in 
vacuum (i.e. T = n = 0). In the language of quantum 
field theory the vacuum model is renormalizable and the 
ultraviolet cutoff scale can be chosen arbitrarily large. 
These arguments show that observables expressed in di- 
mensionless ratios, for example the gap in units of the 
density A/e_p are independent of the choice of A. 



III. SUPERFLUIDITY AND PHASE 
TRANSITION 

The presence of superfluidity is intimately connected 
with the spontaneous breaking of the global continu- 
ous U(l) symmetry of phase rotations. According to 
Goldstone's theorem, a spontaneously broken continu- 
ous symmetry gives rise to the existence of a masslcss 
bosonic mode. The vanishing mass scale m^, — > of this 
mode in turn causes a diverging spatial correlation length 
£ ~ 1/ m v — > oo. This long range coherence then provides 
for typical signatures of superfluidity, such as frictionless 
flow and the existence of vortices as topological defects. 

In our formalism, we use the spontaneous breaking of 
the global U(l) symmetry as an indicator for the onset of 
superfluidity The effective action formalism can directly 
be used to classify the phases of the system according to 
the symmetries of the thermodynamic equilibrium state, 
giving access to a universal description of the phase tran- 
sition. 

In a homogeneous situation, we can consider the effec- 
tive potential U, which corresponds to the average poten- 
tial Uk for k — > or, more precisely for /c _1 equal to the 
macroscopic size of the experimental probe. The effective 
potential can only depend on the invariant combination 
p = ip*ip. The equation of motion then reduces to a 
purely algebraic relation reflecting the condition that the 
effective potential has an extremum at the equilibrium 
value po, 



dU 
d^* 



I'u 



dU 
~dp~ 



Pu 



■ <A> = m l(Po) • ¥>o = 0. 



(26) 



bosonic mass term" rat, as the p 



We have defined a 

derivative of the effective potential. If nonzero, this mass 
term provides a gap in the effective boson propagator, 
and thus suppresses the boson propagation with a spatial 
correlation length £ ~ l/m v . This simple equation can 
be used to classify the phases of the system, 



Symmetric phase : 
Symmetry broken phase : 
Phase transition : 



Pa = 0, 
Pa > 0, 
Pa = 0, 



m; > 0, 



p = 0, 
I = 0, (27) 

where po denotes the solution of Eq. (f2"r?|) . In the sym- 
metric phase (SYM), we deal with a normal gas: there 
is no condensate, po = 0, and the bosonic mass does not 
vanish. The phase with spontaneous symmetry breaking 
(SSB) is instead characterized by a nonvanishing conden- 
sate po > 0. This requires the vanishing of the mass term 



m 2 : Goldstone's theorem is directly implemented in the 
effective action formalism. The massless Goldstone mode 
is responsible for superfluidity, and the vanishing mass 
is associated with a diverging spatial correlation length 
for the Goldstone mode. Without loss of generality we 
choose the phase of spontaneous symmetry breaking by 
fixing fxfi = \/2pa, </?2,o = such that the Goldstone 
mode corresponds to y>2(x). The second "radial mode" 
<pi(x) — <pi,o describes a field with a finite correlation 
length 



in = v/l/(2pot/"(Po)) 



(28) 



(here, we have neglected possible differences in the wave 
function renormalization for the Goldstone and radial 
modes in the SSB regime). 

The flow equation in the SYM regime 



d t m 2 „ = d. 



dUk 
Op 



Po=0 



d tPo = 0, (29) 

has already been discussed before, cf. Eq. (|2"5)l . In the 
SSB phase, the Goldstone mass m 2 = (dU/dp)(po) = 
vanishes identically. The flow equation for the conden- 
sate follows by imposing this minimum condition at any 
scale in the SSB regime, 







d 
It 



dUk 
dp 



= d, 



dU k 



=> dtPo : 
with quartic coupling 



dp 

dtU' k (p ) 



dtPo 



d 2 U k 
dp 2 



K 



K 



U' k '( P o) 



(30) 



(31) 



The phase transition is characterized by the simultaneous 
vanishing of the mass term and the condensate for k — > 0. 
This additional constraint allows to extract the critical 
temperature T c from Eq. (|27j). 

Let us discuss the emergence of spontaneous symmetry 
breaking in our formalism. Independently of the strength 
of the interaction, at low enough temperature the flow of 
the mass term (|29l) hits zero at a finite value of the flow 
parameter k. At this point we switch from Eq. (|29|) to 
Eq. (|30|) . The flow then continues in the SSB regime, 
and a condensate builds up. At fixed c = ak-p and for 
T < T c , the condensate saturates at some positive value 
Po for t — > — oo, i.e., when the cutoff is removed. This 
indicates the presence of spontaneous symmetry break- 
ing. In contrast, for T > T c , the flow drives po back to 
zero at a finite cutoff scale k = fcgri. ln this case, we have 
to switch back to the SYM formulae (|30[) . The presence 
of a nonzero po{k) in the range /csr < k < fcssB can 
be interpreted as the formation of local order on length 
scales between feg SB and fcg R , which is then destroyed by 
fluctuations which persist on larger length scales. 

The scale fcgsB where SSB occurs in the flow varies 
strongly throughout the crossover: It is instructive to dis- 
cuss SSB in both the BEC and the BCS regime. In the 
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BEC region, SSB appears already at a high scale k. The 
condensate builds up quickly during the flow. At very low 
temperature T « 0, it represents the dominant contribu- 
tion to the particle density, apart from a tiny condensate 
depletion which is consistent with a phenomenological 
Bogoliubov theory for bosonic molecules (described by a 
scalar field tp), with a four-boson interaction extracted 
from the solution of the vacuum problem. In the next 
section we will compute the effective molecule scatter- 
ing length om and we find am = 0.72a. This demon- 
strates the relevance of the inclusion of bosonic vacuum 
(or quantum) fluctuations - omitting them would result 
in a Bogoliubov theory with sm = 2a. At higher temper- 
ature, we find that bosonic fluctuations tend to lower the 
value of the condensate at smaller k. We find a second 
order phase transition throughout the whole BEC-BCS 
crossover, corresponding to the expected second order 
phase transition of an 0(2) model. This behavior may 
be inferred from Fig. [4j where the behavior of the gap 
parameter is studied as a function of temperature. One 
clearly observes the continuous closing of the gap as the 
critical temperature is approached. Omitting the bosonic 
fluctuations, there would be no mechanism driving the 
condensate down quickly enough - the phase transition 
then seems to be of first order. This artefact often oc- 
curs in other approaches to the crossover problem at fi- 
nite temperature pjl [HI, [33| and is most severe in the 
BEC regime. This demonstrates the importance of the 
inclusion of effective thermal or statistical bosonic fluc- 
tuations. 

On the BCS side, condensation only appears in the 
deep IR flow for k <C Auv- This goes in line with the 
fact that already very low energy scales like a small tem- 
perature, destroy condensation or pairing in this regime. 
Condensation is indeed a tiny effect in the BCS regime. 
Only an exponentially small fraction of fermions around 
the Fermi surface can contribute to pairing, as seen in 
Fig. H 

Our construction demonstrates the emergence of an ef- 
fective bosonic theory in two respects: The first one con- 
cerns the effective bosonic theory for the tightly bound 
molecules in the BEC regime. The ground state at T = 
is indeed a Bose-Einstein condensate, with a small deple- 
tion due to bosonic self-interactions. Increasing the tem- 
perature, bosonic fluctuations drive the system back into 
the symmetric phase. This second step is precisely the 
mechanism observed previously in purely bosonic theo- 
ries [3 Hll . The presence of bosonic long range fluctu- 
ations is illustrated in Fig. \E\ where we plot the corre- 
lation length which relates to the bosonic mass term as 
indicated in Eq. (|28|) . This aspect of bosonic behavior 
extends over the whole crossover and is discussed in more 
detail in Sect. EH 

In this context, we note a shortcoming of the present 
truncation in the bosonic sector which manifests itself 
in the deep infrared regime of the flow in the low tem- 
perature phase. For example, in the present truncation 
the four-boson coupling X v vanishes polynomially at zero 



A/ £F 



1.5 - 
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FIG. 4: Gap in units of the Fermi energy A/e F = h v -J~po/ 'ep 
as a function of temperature T/Tp. We show the curves ob- 
tained on the BCS side with c _1 = —1 (dashed line), in the 
unitary regime with c _1 = (solid line) and on the BEC 
side with c _1 = 1 (dashed-dotted line) and c _1 = 2 (dotted 
line). As expected the gap for the fermions becomes larger 
on the BEC side for small temperatures. The results clearly 
establish the second order nature of the phase transition. 



temperature, while a more sophisticated infrared analy- 
sis for a weakly interacting Bose gas reveals a logarithmic 
flow of this coupling [36H38j | . In fact, the wave function 
renormalization Z v , i.e. the coefficient to the linear fre- 
quency term in the inverse boson propagator, flows to- 
wards zero. In a more accurate extended truncation a 
quadratic frequency term is generated by the renormal- 
ization group flow. The thereby modified power counting 
results in logarithmic rather than polynomial divergen- 
cies. The correct behavior could thus be implemented 
by extending the truncation with a quadratic frequency 
term in the inverse boson propagator [28|, H3, [38|- We 
note, however, that our results for thermodynamic quan- 
tities are not affected by the deep infrared physics at 
much longer wavelengths than the de Broglie wavelength 
or the mean interparticle spacing. Furthermore, the crit- 
ical behavior at the finite temperature phase transition 
cannot be affected by such an extension of the truncation, 
as only the Matsubara zero mode governs the behavior 
in this region. 



IV. VACUUM LIMIT 

In this section, we consider the effective action in the 
vacuum limit, which is obtained from m the limit 

n — > 0, T — > 0. In this limit, the effective action gener- 
ates directly the basic building blocks for the amputated 
7i-point scattering amplitudes, that in turn directly yield 
the cross section Only simple tree diagrams have 

to be evaluated in order to compute, for example, the 
two-atom scattering via boson exchange. For more com- 
plicated scattering observables where loop effects enter, 
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FIG. 5: Radial correlation length or healing length in units 
of the Fermi momentum ^nkp = ^A^po) - &f as a function 
of temperature T/Tp. We show the curves obtained on the 
BCS side with c -1 = —1 (dashed line), in the unitarity regime 
with c _1 = (solid line) and on the BEC side with c _1 = 1 
(dashed-dotted line) and c _1 = 2 (dotted line). 



jU — em/2 




- 2 2 4 

c" 1 



FIG. 6: Chemical potential for T — minus half the bind- 
ing energy £m/2 = — #(c -1 )c~ 2 as a function of the inverse 
concentration c _1 = (akp)^ 1 . We compare our FRG re- 
sult (solid) to extended mean field theory (dotted; no bosonic 
loops) and a result based on the Schwinger-Dyson equations 
(dot-dashed) [3. 

one observes a considerable simplification of the diagram- 
matic structure of the flow equations in the vacuum limit. 

In order to make contact with experiment, we have 
to relate the microscopic parameters which characterize 
the theory at a high momentum scale Atjv to the macro- 
scopic observables for two-atom scattering in vacuum, 
like the scattering length a or the molecular binding en- 
ergy. This is the analogue of renormalization in quantum 
electrodynamics, which relates the bare (microscopic) to 
the renormalized (macroscopic) coupling. For our set- 
ting the bare couplings diverge in the limit of infinite UV 
cutoff Atjv, while all observable quantities remain finite 
once they are expressed in terms of the renormalized pa- 



rameters. In our approach, this mapping between bare 
and renormalized parameters is implemented by choosing 
the initial parameters at the microscopic scale Atjv such 
that a few two-body observables in vacuum are matched 
to their experimental values in the limit k — > 0. 

In our picture, atom scattering in vacuum is mediated 
by the formation and dissociation of a collective boson. 
We stress, however, that in the limit of broad resonances, 
h v — > oo, the collective bosonic degree of freedom is 
purely auxiliary on the level of the microscopic action 
(defined at Atjv) [3- In this case the microscopic in- 
teraction is simply a pointlike contact interaction. Nev- 
ertheless, fcrmionic fluctuations dynamically generate a 
physical bosonic bound state on the BEC side of the res- 
onance where the fermionic scattering length is positive. 
For k — > 0, we find bosonic degrees of freedom with dis- 
persion uj = q 2 /(AM) - the boson mass 2M is appropri- 
ate for the composite objects. On the BCS side close to 
the Feshbach resonance (negative a) the bosonic particle 
remains as a resonant state, with energy now somewhat 
above the open channel energy level. 

In the broad resonance limit and not too far away 
from the Feshbach resonance, the physics of both the vac- 
uum and the thermodynamics becomes very insensitive 
to details of the microscopic parameters [To|, [HI [2l[ . We 
find that further microscopic information is suppressed at 
least ~ h~ 2 for large h v . In the RG language, this "uni- 
versality" is conveniently expressed in terms of a (non- 
trivial) broad resonance fixed point. As a signature of 
this strong universality, the two-channel model with ex- 
plicit molecular degrees of freedom becomes equivalent 
to a single channel model, in which the interaction is a 
priori specified by a single scale, the scattering length - 
both models are characterized by the same fixed point. 
We recover this feature directly from the solution of the 
flow equations. 

The physics of the vacuum limit describes few-body 
scattering, which can also be computed from quantum 
mechanics. The functional integral approach is equiva- 
lent to quantum mechanics, and the flow equations are 
exact. The results from a solution of the flow equations 
for k — > should therefore coincide with the quantum 
mechanical results, up to artifacts of the truncation. One 
may employ the comparison to quantum mechanics for 
an estimate of the error due to the truncation. 

At first sight, the use of the flow equation machinery 
for the solution of a quantum mechanical problem may 
seem unduely complicated. We recall, however, that the 
flow equation result can immediately be extended to the 
many-body problem at nonzero density or temperature, 
simply by changing the parameters jl and T in the flow 
equation (|16j) . This is not possible for quantum mechan- 
ics. Furthermore, many different microscopic Hamilto- 
nians can now be investigated with comparatively little 
effort - it suffices to change the initial values of the flow. 
The fixed-point behavior offers an easy access to uni- 
versality, which cannot be seen so directly in quantum 
mechanics. 
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The vacuum problem can be structured in the follow- 
ing way: The two-body sector describes the pointlike 
fcrmionic two-body interactions. It involves couplings up 
to fourth order in the fermion field in a purely fermionic 
setting. In the language with a composite boson field 
if ~ ipif) we have in addition terms quadratic in ip or 
~ ipipip* + h.c. The two-body sector decouples from 
the sectors involving a higher number of particles, de- 
scribed by higher-order interactions parameters such as 
the dimer-dimer scattering (which corresponds to inter- 
actions of eighth order in tp in a fermionic language) . This 
decoupling reflects the situation in quantum mechanics, 
where a two-body calculation (in vacuum) never needs 
input from states with more than two particles. 

Not too far from unitarity the solution of the two-body 
problem fixes the independent renormalized couplings 
completely. The couplings in the sectors with higher 
particle number then are derived quantities which can 
be computed as functions of the parameters of the two- 
body problem. Here, we will study the dimer-dimer or 
molecular scattering length am as an important example. 

We emphasize that this almost complete "loss of mi- 
croscopic memory" near a Feshbach resonance is char- 
acteristic for our system with two species of fermions. 
For three species with SU(3) symmetry, or for bosonic 
Feshbach resonances, the three-body sector will involve 
an additional "renormalized coupling" , associated to the 
energy of a trion bound state in the SU(3) case (39l. |40|. 
For the present system, however, we expect that close to 
a (broad) Feshbach resonance all "macroscopic" physi- 
cal phenomena can be expressed by only one renormal- 
ized dimensionlcss parameter, namely the concentration 
c = akp, kp = [i^n) 1 ^ , plus the dimesionless temper- 
ature. The inverse concentration c _1 is a relevant pa- 
rameter in the sense of critical phenomena, with c^ 1 =0 
being the location of the fixed point. Scales may be set 
by the density n or the scattering length a. 



A. Projection onto the Vacuum 

First we specify the prescription which projects the 
effective action onto the vacuum limit [l6| : 



T(vac) 



lim r fc=0 

T,fc F ->0 T>T C 



(32) 



Here kp can be viewed as the inverse mean interparti- 
cle spacing kp ~ 1/d. Taking the limit kp — > then 
corresponds to a diluting procedure where the density of 
the system becomes arbitrarily low. However, the limit 
is constrained by keeping the dimensionlcss temperature 
T = T/e F = 2MTjk\ above criticality. This ensures 
that many-body effects such as condensation phenomena 
are absent. In this limit, our functional integral approach 
becomes equivalent to quantum mechanics. 

We find that for n = T = the crossover at finite 
density turns into a second order "phase transition" in 
vacuum fiol . [l6| as a function of fh 2 (A) which is given, in 



turn, by the magnetic field B. This is expressed by the 
following constraints which separate the two qualitatively 
different branches of the physical vacuum, (cf. Eq. (|27p) 



m 2 \ k - 
mjk- 



>o > 0, /i = atom phase (a -1 < 0) 

>o = 0, /i < molecule phase (a -1 > 0) . (33) 



0, \i — resonance 



These formulae have a simple interpretation: On the BCS 
side, the bosons experience a gap m 2 > and the low- 
density limit describes only fermionic atoms. On the 
BEC side, the situation is reversed: the fermion prop- 
agation is suppressed by a gap —fj,. The ground state 
is a stable molecule, and the fermionic chemical poten- 
tial can be interpreted as half the binding energy of a 
molecule, cm = 2/i [l6| - this is the amount of energy 
that must be given to a molecule to reach the fermionic 
scattering threshold. 

Evaluating the bosonic mass term fh 2 , in vacuum on 
the BEC side (see below), one finds the well-known 
universal relation between binding energy and scatter- 
ing length in vacuum, em = — l/(Ma 2 ) in the broad- 
resonance limit h v — > oo. Close to the resonance, 
em 0, this holds for arbitrary h v as well as other mi- 
croscopic details and hence establishes the second order 
nature of the vacuum phase transition. 



B. Diagrammatic simplifications 

On the technical side, the procedure specified above 
leads to a massive simplification of the diagrammatic 
structure as compared to the finite density and tempera- 
ture system. With the aid of the residue theorem, we can 
prove the following statement: All diagrams with inner 
lines pointing into the same direction ( thereby forming a 
closed tour of particle number) do not contribute to the 
flow in vacuum. Here, the direction of a line is given by 
the go flow of the propagator. 

For a proof, we first consider the form of the renor- 
malized inverse fermion and boson propagators in the 
vacuum limit which there become diagonal and are rep- 
resented by the entries 



Pp(q) = iq Q + q 2 - fx, 
P v (q) = lS v q + q 2 /2 + 



(34) 



Lines pointing into the same direction represent integrals 
over products of PpT^P" 1 with the same sign of the 
frequency variable in the loop. Without loss of generality, 
we can choose it to be positive. 

In the presence of a nonzero cutoff, the spacelikc real 
part of the regularized inverse propagators Pp + , P v + 
i?^, including the mass terms, is always positive. Hence, 
the poles all lie in the upper half of the complex plane. 
Closing the integration contour in the lower half-plane, 
no residues are picked up, implying that these integrals 
vanish. A derivative with respect to the cutoff Rk in- 
crements the number of inner lines Pp 1 , P^ 1 by one (it 
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changes the multiplicity of the poles), but does not affect 
the sign of the real part of the propagator nor of go > such 
that our argument remains valid for both the regularized 
loops and their k derivative entering the flow equation. 

Let us now discuss different types of diagrams explic- 
itly which do not contribute to the vacuum flow. 

• The mixed diagram with one inner fermion and one 
inner boson line, driving the renormalization of the 
fermion propagator, does not contribute. This im- 
plies for all k 

P F (q) = iLU + q 2 -n. (35) 

• The box diagram with two inner boson and fermion 
lines, in principle generating a four-fcrmion inter- 
action even for vanishing initial coupling, does not 
appear. For A^|fc=A = we thus have 

= (36) 

at all scales. This means that partial bosonization 
is very efficient in this limit - fluctuations are com- 
pletely absorbed into the bosonic sector, and there 
is no "backreaction" on the fermion propagator. 

• Contributions to the flow from diagrams involving 
solely one inner line are zero. This is directly re- 
lated to the fact that the particle density vanishes 
by construction - diagrammatically, the trace over 
the full propagator is represented as a closed loop 
coupled to an external current, the chemical po- 
tential (tadpole graph). This is relevant for the 
bosonic contribution to the boson mass which thus 
vanishes - the renormalization of the boson mass in 
the physical vacuum is purely driven by the fermion 
loop. 

• The bubble diagrams (involving two inner boson 
or two inner fermion lines forming a closed tour) 
vanish. 



The ladder diagrams instead, which have lines point- 
ing in the same direction, contribute. Due to the opposite 
signs of the momentum variables in the propagators, the 
poles of these diagrams are located in both the upper 
and the lower half plane. It is worth noting that the 
four-boson coupling receives such a bosonic ladder con- 
tribution (cf. Fig. [3]), which will be discussed in Sect. 
HVGl 

C. Solution of the two-body problem 



The two-body problem is best solved in terms of the 
bare couplings. Their flow equations read 

h 2 

d k Z v = -_|^(fc 2 +M )(fc 2 + M )3/ 2 , 

d k A v = - V A v A v /k = -^0(k 2 + M )(fc 2 + M ) 3 / 2 , 
d k h v = 0. (37) 

The flow in the two-body sector is driven by fermionic di- 
agrams only. There is no renormalization of the Feshbach 
coupling due to U(l) symmetry. This statement holds in 
the absence of an initial coupling, X^,\k=A = 0, where the 
whole scattering amplitude of the fermions is expressed 
by the exchange of a bosonic bound state. (Taking an 
additional background coupling into account leads to a 
renormalization of h v [2l[ which is compatible with the 
U(l) invariance.) 

Eqs. (|37[) are solved by direct integration with the 
result 
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ml{k) 



m 2 (A) - 6(A 2 + ^ 



+0(fc 2 + M )^ 



Z v (k) = Z V (A) - 0(A 2 + [VA^ (5A2 A t 2/i) 



48tt 2 



3 / v/A^TT^ 



arctan 



/T, ( 5fc2 + 2 A0 



3 / vW+Ji 



k 4 



arctan 



-/' 



(38) 



(39) 



A v {k) = A !p (A) - fl(A 2 + M ) [7A^ (5A2 A t ^ 



48tt 2 



J= arctan 



/i 2 r c 5A ,2 + 2 \ 



48tt 2 



A- 4 



j — arctan ( — = ^ 



(40) 



Here, A is the initial ultraviolet scale. Let us discuss the 
initial value for the boson mass. It is given by 



fh 2 {K) = v(B)-2n + 5v(A). 



(41) 



It features a physical part, the detuning v(B) = hm(B — 
Bo), which describes the energy level of the microscopic 
state represented by the field ip with respect to the 
fcrmionic state ip. At a Feshbach resonance, this energy 
shift can be tuned by the magnetic field B, /j,m denotes 
the magnetic moment of the field ip, and Bq is the reso- 
nance position. Physical observables such as the scatter- 
ing length and the binding energy are obtained from the 
full effective action and are therefore related to the cou- 
pling constants at the infrared scale k = 0. The quantity 
5is(A) denotes the renormalization counter term that has 
to be adjusted such that the condition (|33|) is satisfied 
for k — > on the resonance, see below. Equations similar 
to Eq. (|37p are derived in [4l|, |42j], where other cutoff 
functions are used, however. 



D. Renormalization 

We next show that close to a Feshbach resonance the 
"microscopic" parameters (initial conditions to the flow 
equations) m^A = fn^ik = A) and A = /i 2 (fc = A) 
which enter our calculations are related to the physical 
observables, the detuning from the Feshbach resonance 
v = (Im(B — Bo) and its width /i 2 . Alternatively to 
the set {v, h 2 ,}, we may also choose {a,h 2 }, where a is 
the scattering length. The idea behind the procedure of 
fixing the parameters in our RG framework is illustrated 
and discussed in Fig. [TJ The relations are 



and the physical observables are related by 



K,A 



8n m (B - B )' 



(43) 



In the vacuum limit specified above, and in the absence 
of background interactions, the Yukawa or Feshbach cou- 
pling h Vt \ is not renormalized, and so does not change 
with scale. Therefore, unlike the mass term which re- 
ceives UV renormalization, the Feshbach coupling is a 
direct physical observable, associated to the width of the 
resonance (see below). In the limit of a broad Feshbach 
resonance only a single parameter is needed to describe 
the scattering physics - only the ratio of the Feshbach 
coupling /i 2 A and the detuning is physical. Away from 
the Feshbach resonance the Yukawa coupling may de- 
pend on B, h 2 v {B) = /i 2 + d(B - B ) + . . . . Also the 
microscopic difference of energy levels between the open 
and closed channel may show corrections to the linear in- 
dependence, v{B) = fiM(B — -Bo) + c%(B — Bo] 2 + ... or 
Mm -> Mm + c 2 (B - B ) + .... Using h 2 v {B) and fi M (B) 
our formalism can easily be adapted to a more general ex- 
perimental situation away from the Feshbach resonance. 
The relations in Eqs. (|4"2j) and (|4"3")l hold for all chemical 
potentials \i and temperatures T. For a different choice 
of the cutoff function the coefficient 5v(A) being the term 
linear in A in Eq. (|42[) might be modified. 

We want to connect the bare parameters m 2 A and 
n 2 A with the magnetic field B and the scattering length 
a for fcrmionic atoms as renormalized parameters. In our 
units, a is related to the effective interaction A^. c ff by 



A 



Vi.cff 

"8tT' 



(44) 



Hm(B - -B ) - 2[i 



h 2 
6tt 2 



A, h 



(p,A 



In the absence of a background interaction, the fermion 

(42) 

interaction A^,. ff is determined by the molecule exchange 
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process in the limit of vanishing spatial momentum 



V,eff : 



P ip (lu,p 2 = Q,^) 



(45) 



Here, we work with Minkowski frequencies lo related to 
the Euclidean ones by iqo — > — Even though (|45|) 
is a tree-level process, it is not an approximation, since 
P v = P v \k->o denotes the full bosonic propagator which 
includes all fluctuation effects. The frequency in Eq. (|4"5)) 
is the sum of the frequency of the incoming fermions 
which in turn is determined from the on-shell condition 



lu = 2lu^ = — 2(1. 



(46) 



On the BCS side we have n 
find with 



(see Eq. (J33J) and 



P v (uu = Q,q = Q) =m 2 v {k 



0) 



(47) 



the relation 



On the BEC side, we encounter fi < and thus lu > 0. 
We therefore need the bosonic propagator for lu ^ 0. 
Even though we have computed directly only quantities 
related to P v at lu = and derivatives with respect to lu 
(Z^, we can obtain information about the boson propa- 
gator for nonvanishing frequency by using the scmilocal 
U(l) invariance described in App. [B] In momentum 
space, this symmetry transformation results in a shift of 
energy levels 



-> <p(lu-26,p) 
fx — > /i + <5. 



(53) 



Since the effective action is invariant under this symme- 
try, it follows for the bosonic propagator that 



P v (u,p,p) = P<p(lu - 25, p, p + 5). 



(54) 



To obtain the propagator needed in Eq. (|4"5j) , we can use 
5 = —fi and find as in Eq. (|52l) 



V,eff 



ft 2 



(48) 



where m 2 = m^(i = 0). For the bosonic mass terms 
at fi = 0, we can read off from Eqs. (|40|) and (|41|) that 



-*£A - M m(B - So) + ^(A) - -P^-A. 



6tt 2 
(49) 

To fulfill the resonance condition in Eq. (|33|) for B = Bq, 

p = 0, we choose 



(Si/(A) 



6tt 2 



A. 



(50) 



The shift 5v (A) provides for the additive UV rcnormaliza- 
tion of to 2 as a relevant coupling. It is exactly cancelled 
by the fluctuation contributions to the flow of the mass. 
This yields the general relation (|42|) (valid for all p) be- 
tween the bare mass term m 2 A and the magnetic field. 



On the BCS side we find the simple vacuum relation 



"4,0 = Mm(-B - Bq). 



(51) 



Furthermore, we obtain for the fermionic scattering 
length 



8n m (B-B )' 



(52) 



This equation establishes Eq. (|43|) and shows that ft 2 A 
determines the width of the resonance. We have thereby 
fixed all parameters of our model and can express m 2 
and ft 2 A by B 



B Q and a. The relations (O and (|35)l 
remain valid also at nonzero density and temperature. 
They fix the "initial values" of the flow (ft 2 — > ft, 2 A ) 
at the microscopic scale A in terms of experimentally 
accessible quantities, namely B — Bq and a. 



A 



ft 2 



ft 2 



P v (cu = 0,jP =0,/x = 0) 



Mm(-B — ^o) 

(55) 

Thus the relations (|42p and (f4"5)l for the initial values 
m^A and ft 2 A in terms of B — Bq and a hold for both 
the BEC and the BCS side of the crossover. 



E. Binding energy 

We next establish the relation between the molecu- 
lar binding energy em, the scattering length a, and the 
Yukawa coupling ft, 2 A . From Eq. (|4U|) . we obtain for 



k = and n < 



^3,0 



/i M (-B - -B ) - 2p 
ft 2 

V,A 



6tt 2 



A- ^ + 5(1-^) 



In the limit A/, 



i^arctan ^ 



co this yields 



//m(-B - -B ) - 2/x 



Together with Eq. (|52]l. we can deduce 

ft 2 



87 r« + 2M-^^ 



8tt 



(56) 



(57) 



(58) 



which holds in the vacuum for all p. On the BEC side 



where 



this yields 



(59) 
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The binding energy of the bosons is given by the dif- 
ference between the energy for a boson m 2 / 'Z v and the 
energy for two fcrmions 

l tp,a 



-2fi. On the BEC side, we can 



2/x 



(60) 



From Eqs. ([55]) and (f6"0"|) we find a relation between the 
scattering length a and the binding energy cm 



CM J 



i (8tt) 2 



(61) 



our findings can easily be taken over to this more complex 
situation. 

Let us consider the flow equations (|3T[) for the renor- 



malized quantities (m 2 = rn^/A^ 
Ztp/Atp) in the regime where k 2 ^> —fj, 



h 2 = 



K/a 



Sip — 



dt < = 6^ 



d t S v 



6ir 2 k 
d t A v 



——h 2 

6ir 2 k v ' 



(63) 



In the broad resonance limit /i 2 A — > oo, this is just the 
well-known relation between the scattering length a and 
the binding energy cm of a dimcr 



1 



£M = 9 

a 1 



Ma 2 ' 



(62) 



The last two terms in Eq. (|61[) give corrections to Eq. 
(|62[) for more narrow resonances. Only in this case, the 
coupling ^ acquires an independent physical meaning; 
for broad resonances, the only physical observable is the 
scattering length a in Eq. (|43|) . 

The solution of the two-body problem turns out to be 
exact as expected. In our formalism, this is reflected 
by the fact that the two-body sector decouples from the 
flow equations of the higher-order vertices: no higher- 
order couplings such as A^, enter the set of equations 
([37)) . Inspection of the full diagrammatic structure of 
our truncation shows that there is a tadpole contribution 
~ X v at arbitrary n and T which enters the flow equation 
for the bosonic mass. However, the contribution of the 
corresponding diagram to the flow vanishes in vacuum. 
Extending the truncation to even higher order vertices 
or by including a boson- fermion vertex i\vi\)ip*ip does not 
change the situation, since there are only two external 
lines in the two-body problem and the flow equation in- 
volves only one-loop diagrams, such that contributions 
from these vertices do not appear. We emphasize that 
this argument is not bound to a certain truncation. 



F. Universality 

Universality means that the macroscopic physics (on 
the length scales of the order of the inter-particle dis- 
tance) becomes independent of the details of the micro- 
physics (on the molecular scales) to a large extent. This 
is due to the presence of fixed points in the renormaliza- 
tion flow. Approaching the fixed point, the flow "loses 
memory" of most of the microphysics, except for a few 
relevant parameters. We will show here how the broad 
resonance universality emerges within the flow equations 
for the vacuum. The fixed point structure remains simi- 
lar for non-vanishing density and temperature, such that 



In this scaling regime, the flow loses memory of all initial 
conditions (except for the mass term which is a relevant 
parameter) , provided that the dimcnsionless combination 
h^/h is attracted to a non-perturbative fixed point. The 
flow equation 




hi 



(64) 



exhibits indeed an IR-attractive fixed point (scaling solu- 
tion) given by t]a. v = 1, h/^/k = 6ir 2 . This fixed point is 
approached rapidly provided the initial value of ft 2 (A)/ A 
is large enough. Memory of the precise value of ft. 2 , (A)/ A 
is then lost. The situation is similar for most other pos- 
sible additional couplings that we have not displayed 
here: The dimcnsionless renormalized couplings are all 
attracted to fixed point values. For example, this also 
applies to the four-boson coupling X v [2lj |. 

A notable exception is the mass term to 2 . The dimcn- 
sionless combination obeys 



1,1 I 67T 2 k ^ ^ k 2 k 2 ' 



where the fixed-point value for /i 2 /fc has been used in 
the second equation. The fixed point at m 2 /fc 2 = 1 is 
unstable for the flow towards the infrared. A small ini- 
tial deviation from the fixed-point value tends to grow 
as k is lowered. The mass term is therefore a relevant 
parameter. For k — > 0, to 2 approaches a finite value, see 




Eq. (|38|) , and sets the scale for all dimensionful quantities 
on the BCS-side of the crossover. 

We conclude that the precise initial values of all quan- 
tities except for the mass term are unimportant. For 
definitcness, we choose A v = Z v = 1, X v = 0. 

Following the k evolution further, on the BEC side 
for fi < the flow leaves the scaling regime near k ~ 
y/—fi. Thus a negative /i also plays the role of a relevant 
parameter, similar to m 2 for /i = 0. The IR limit k = 
is determined by a single scale set by the value of jj,. This 
reflects the single relevant parameter which characterizes 
the infrared physics of a broad resonance on the BEC 
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side. Wc consider the flow of S v with h^/k assuming its 
fixed point value, 



l )VA v - 



(65) 



This equation is solved by S v = Z v /A v = 1 for k — > 0, 
corresponding to the dispersion relation for particles of 
mass 2M . 

To summarize, the flow starts at some microscopic 
scale A with a fine-tuned choice of to 2 A > given by 



Eq. fl42 



For broad resonances, the precise initial values 
of all other parameters of our truncation (S v , X v , . . . ) 
are unimportant as a consequence of universality - the 
only relevant parameter is the inverse scattering length. 
In particular we see how a molecular bound state ap- 
pears on the BEC side of the crossover, independently of 
the precise form of the microscopic physics which may be 
described by a simple pointlikc interaction. 



G. Dimer-Dimer Scattering 

So far we have considered the sector of the theory 
up to order tp*ipil), which is equivalent to the fermionic 
two-body problem with pointlike interaction in the limit 
of broad resonances. Higher-order couplings, in par- 
ticular the four-boson coupling A^,^*^) 2 , do not cou- 
ple to the two-body sector. Nevertheless, a four-boson 
coupling emerges dynamically from the renormalization 
group flow. In vacuum we have po = and X v is defined 
as X v = U%(0), cf. Eq. (HI])- The flow equation for X v 



can be found by taking the p-derivative of Eq. ([22 
V2k 3 ( 2 



kdkX^ = 2r\ A U'l - 



x2(C/0 2 



1 - 



'B,Q 



3s 



(0,1) 
B,Q 



-,va v 



h V J 2 ) 
7& I 



V2k 5 X 2 v 



3tt 2 S v (mj + fc 2 ) 2 



h%e{p + k 2 ){p + k 2 ) 3 ' 2 



3ir 2 k 3 F '' 
(1-2 VA J5) 

(66) 



47T 2 fc 6 

There are contributions from fermionic and bosonic vac- 
uum fluctuations (cf. Fig. [3]), but no contribution from 
higher p derivatives of U as can be checked by apply- 
ing the arguments from Sect. IIV Bl The fermionic dia- 
gram generates a four-boson coupling even for zero initial 
value. This coupling then feeds back into the flow equa- 
tion via the bosonic diagram. 

First, we discuss the implications of broad-resonance 
universality for the four-boson coupling and the ratio of 
the scattering length am for molecules and a for fermions. 
In the scaling regime k 2 ^> — p and at large ft. 2 , we use 
the scaling form Eq. (|63[) implying h\ 



v - 6n z k, j] v 



1, m 2 = k 2 . The dimensionless ratio Q = X v k 3 /h 

9\/2V ^ 



obeys the flow equation 
kd k Q = 3Q 



1 

i^ 2 



(67) 



The flow exhibits an infrared stable fixed point. Its value 
Q* ~ 0.008 corresponds to a rcnormalizcd coupling scal- 
ing ~ k~ x , 



367r 4 Q >t 



(68) 



This can be compared to the effective four-fermion cou- 
pling in the scaling regime, A^, e ff = —h^/m 2 , = —6w 2 /k. 
Thus A^/A^eff approaches a fixed point. The constant 
ratio between these two quantities in the scaling regime 
of the flow is the origin of the universal ratio between the 
scattering length for molecules and atoms. 

Having settled the universality of the scattering length 
ratio, we now address its numerical value. The scattering 
lengths are related to the corresponding couplings by the 
relation (cf. [H|) 



9 

— \ ' 



A 



V>,eff 



8ttc 



8tt 



(69) 



The last equality is obtained in the molecule phase for 
p < and k = 0, using Eq. ([56]) together with the con- 
straint (|33[) . Omitting the molecule fluctuations, a direct 
integration of Eq. (|66[) yields A v = 8tt/^/— p and there- 
fore ajvi/a = 2. Inclusion of the molecule fluctuations 
lowers this ratio. We stress that the vacuum value for 
A^/A^eff has to be computed for k — > 0. It therefore dif- 
fers somewhat from the ratio in the scaling regime (Eq. 
([55)) ) since the final flow for k 2 < —p has to be taken into 
account. The result is weakly cutoff dependent. With our 
truncation and choice of cutoff one finds ayi/a = 0.718. 

The ratio a^i/a has been computed by other meth- 
ods. Diagrammatic approaches give om/s = 0.75(4) [43| . 
whereas the solution of the 4-body Schrodinger equation 
yields om/i = 0.6 [44|, confirmed in QMC simulations 
[45l | and with diagrammatic techniques [46| . Our calcu- 
lation can be improved by extending the truncation to 
include a boson-fermion vertex X v ^ which describes the 
scattering of a dimer off a fermion [4l|. Inspection of 
the diagrammatic structure shows that this vertex in- 
deed couples into the flow equation for X v . Moreover, 
this coupling is important for a precision estimate of the 
dynamically generated atom-dimer scattering length. 



V. PARTICLE-HOLE FLUCTUATIONS 

In the previous section we were concerned with the 
flow equations in the vacuum limit where both the den- 
sity and the temperature vanish, n = T = 0. It is one 
of the main advantages of our method that we can treat 
this limit as well as the many-body problem in thermal 
equilibrium with nonzero density and temperature in a 
unified approach. On large scales k 2 ^> T, k 2 ^> n 2 / 3 the 
flow equations are basically the same as in the vacuum 
case. The qualitative features and the fixed points are 
the same as discussed in the previous section. However, 
for scales of the order of the inverse particle distance 
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k sa n 1 / 3 the flow equations arc modified by many-body 
effects. We have seen in Sec. Illll how the fermionic fluctu- 
ations lead to a nonvanishing order parameter po > for 
k < fcssB connected to local order, and to superfluidity 
if pa remains positive for k = 0. 

Other interesting many-body effects are the screen- 
ing of the four-boson interaction X v on the BEC side 
of the crossover or the effect of particle-hole fluctuations 
on several thermodynamic observables on the BCS side. 
The first effect is well captured in our present truncation. 
It manifests itself by additional terms in the flow equa- 
tion for Xtp which are not present in the vacuum limit. 
These bosonic fluctuations lead to a decrease of the four- 
boson interaction \ v {k = 0). This is very similar to 
the screening effect in a system of non-relativistic bosons 
with pointlike interaction [28|, H?} . The effect of particle- 
hole fluctuations is more involved. The truncation dis- 
cussed in Sect. Ill Cl is not yet sufficient to capture it. We 
next discuss the extensions necessary to include this ef- 
fect following Ref. [23[ . Within the approach discussed in 
Sect. HIl the effect of particle-hole fluctuations generates 
a four fermion vertex A^, for k < A, even if we start with 
A^a = due to bosonization of the original pointlikc in- 
teraction. The total effective four-fermion interaction is 
composed of a contribution from the exchange of a boson 
and of generated by the renormalization flow from the 
box-diagrams involving the fermions and the bosons as 
displayed in Fig. [3 



leads to a modified flow equation for the Yukawa coupling 
h v . For details we refer to Ref. [23j . 

Most prominently, the effect of particle-hole fluctua- 
tions is seen in the critical temperature T c . For non- 
relativistic fermions the task of a precise computation of 
T c is difficult, the basic reason being that the phase tran- 
sition itself is a non-pcrturbative phenomenon, linked 
to an effective four-fermion interaction A^, je fi growing to 
large values. Indeed, the Thouless criterion states that 
a phase transition is connected with a diverging effective 
interaction between fermions A^, ;C ff — > oo. In order to lo- 
cate T c precisely, one has to follow the scale-dependence 
of A^gff with sufficient precision. The problem arises 
from the substantial momentum dependence of the fluc- 
tuation contributions to the effective four-fermion ver- 
tex. At weak coupling, the dominant deviation from the 
BCS result for T c - which is obtained by neglecting the 
momentum dependence of the four-fermion vertex com- 
pletely - comes from the in-medium screening due to the 
particlc-holc fluctuations. It was first shown by Gorkov 
and Melik-Barkhudarov [4^] (see also (50[) in a pertur- 
bative setting that this effect leads to a decrease of the 
critical temperature in the BCS regime by a factor ~ 2.2 
compared to the result obtained in the original BCS the- 
ory. 

We show the numerical result of Ref. [23| for T c /Tp in 
Fig. [51 together with the critical temperature obtained 
in the truncation (TlH). The results from the extended 



A^off — — -pf- + \ 

Flo 



(70) 



The vertex A^ is generated only in medium, i.e. at 




FIG. 7: Box diagram for the flow of the four-fermion interac- 
tion. Arrows denote the flow of the particle number. 

finite temperature and density while the corresponding 
diagram vanishes in the vacuum. 

If we contract the dashed boson exchange lines in Fig. 
13 to a point it is easy to see that it describes the in- 
teraction of two incoming particles (with momenta pi, 
P2 and similar for the outgoing channel) with a (virtual) 
particle-hole pair. Using the formalism of rcbosoniza- 
tion developed in [48| it is possible to adjust the boson 
exchange in a scale-dependent manner such that the ef- 
fect from particle-hole fluctuations is included in the bo- 
son exchange description. This method involves scale- 
dependent fields where the boson is an effective degree 
of freedom. Technically, the use of scale-dependent fields 
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FIG. 8: (Color online) Dimensionless critical temperature 
T c /Tf as a function of the inverse concentration c _1 = 
(o/cf) 1 . This plot is taken from Ref. [2^]. The black solid 
line includes the effect of particle-hole fluctuations. The result 
obtained when particle-hole fluctuations are neglected is also 
shown (dot-dashed line). For comparison, we plot the BCS 
result without (left dotted line) and with Gorkov's correction 
(left dashed). On the BEC side with c _1 > 1 we show the 
critical temperature for a gas of free bosonic molecules (hori- 
zontal dashed line) and a fit to the shift in T c for interacting 
bosons, AT C ~ c (dotted line on the right). 

truncation used in [23| are depicted be the solid line in 
Fig. El This result agrees with the prediction by Gorkov 
and Melik-Barkhudarov in the regime with small nega- 
tive scattering length (long-dashed line in Fig. [5]). By 
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contrast, the simpler truncation described in Sec. [TT1 of 
this paper (dot-dashed line) reproduces the original BCS 
result (dotted line in Fig. [8]). Both our approximations 
yield the same result for c _1 > 0.5 where the effect of 
particle-hole fluctuations disappears. This is expected 
since the chemical potential is negative here, [i < 0, and 
there is no Fermi surface any more. In the BEC limit 
for very small positive scattering length we find that our 
result approaches the critical temperature of a free Bose 
gas 



T, 



c,BEC 



0.218. 



(71) 



For c — > 0+ this value is approached in the form 



f c ,BEC 



T, 



c,BEC 



1/3 a M 
KQM n M = K 



a (6TT 2 ) 1 / 3 ' 



(72) 



Here, tim = n/2 is the density of molecules and «m is the 
scattering length between them. For the ratio ajj/a we 
use our result ay^/a = 0.718 obtained from solving the 
flow equations in vacuum, see sect. IIV[ since for tightly 
bound molecules this is the relevant microscopic interac- 
tion parameter. For the coefficients determining the shift 
in T c compared to the free Bose gas we find k = 1.55. 
Given the simplicity of our approach, it is remarkable 
that our result is near the value n ~ 1.3 found in an ef- 
fective three dimensional "classical" bosonic theory [fnj 
with lattice simulations [52l |53| and with the functional 
RG using elaborate momentum-dependent truncations 
A direct quantitative comparison is, however, inhib- 
ited by the fact that our result may receive contributions 
from the residual fermionic modes. 

At the unitarity point c _1 = we find T C /T F = 0.264 
with the approximation used in |23| . The simpler trunca- 
tion used in Sect [TT] gives T C /T F = 0.276. This should be 
compared with results from Quantum Monte Carlo sim- 
ulations: T C /T F = 0.15 and T C /T F = 0.245 @. The 
measurement by L. Luo et al. [55j in an optical trap gives 
T C /T F = 0.29(+0.03/ - 0.02), which is a result based on 
the study of the specific heat of the system. Further the- 
oretical predictions based on a variety of methods span 
a whole range of values: T C /T F = 0.249 (e expansion 
0), T C /T F = 0.183 (Borel-Pade approximation [ill), 
T C /T F = 0.136 (1/N expansion [l3|), T C /T F = 0.160 
(2PI techniques [ill]). This demonstrates that quantita- 
tive precision requires a control of the systematic errors 
in a given approximation scheme. A particular important 
point concerns the quantitative accuracy in the determi- 
nation of the density. We will discuss this issue in Sect. 
IVIII If we replace our determination of the density by 
the "free particle density" we obtain (in the truncation 
of Sect. HJ) T C /T F = 0.171. 



VI. UNIVERSAL LONG-RANGE PHYSICS 
CLOSE TO THE PHASE TRANSITION 

As the phase transition is approached, one expects to 
find universal critical behaviour if the temperature is 
close enough to the critical temperature T c . This fea- 
ture can be seen directly from the nonperturbative flow 
equations. The universal behaviour arises from long- 
range bosonic fluctuations with momenta below T 1 / 2 and 
larger than the inverse correlation length £ , and re- 

1 /2 

quires £T C ' ^> 1 . The temperature acts as an effective 
infrared cutoff for the influence of the fermion fluctu- 
ations on the long-range physics. Correspondingly, for 
fc/V^c *C 1 the fermion contribution to the flow becomes 
small and can be neglected. This can be inferred directly 
from the dependence of the fermionic threshold functions 
on k/y/T, as we have discussed after Eq. (TT9")) . On the 
other hand, for k/VT <C 1 the bosonic threshold func- 
tions are transmuted to the ones of a classical statistical 
system for bosons with U(l) symmetry. The infrared flow 
is governed by a Wilson-Fisher fixed point for the N = 2 
universality class, where N = 2 corresponds to the two 
real scalar fields on which global S'0(2)-transformations 
act linearly. Following the flow equations to k — > leaves 
no doubt that this universality class governs the critical 
exponents and amplitude ratios for T sufficiently close to 
T c . For more elaborate truncations, the non-perturbative 
flow eq uations yield precise values for these quantities 

sain. 

It is not our aim here to reproduce the precise univer- 
sal quantities by an extension of our truncation. We are 
rather interested in another question: how extended is 
the temperature range for the universal behavior? Fur- 
thermore, we want to compute the non-universal aspects 
of the phase transition as well. They will depend, of 
course, on the precise location of the system in the BCS- 
BEC crossover. As an example, we investigate the tem- 
perature dependence of the boson self- interaction X v . 
Close to the critical temperature, X v is known to van- 
ish (A^/m 2 goes to a constant in the symmetric phase). 
Since also the gap A vanishes for T — > T c , we may take 
the size of the gap A(T) as a measure of the distance 
from the critical temperature. In Fig. |H]we plot the ratio 
\ V {T)/\ V {T = 0) as a function of A(T)/A(T = 0). The 
plots correspond to the BCS (solid; c _1 = —2), crossover 
(dashed dotted; cT 1 = 0) and BEC (dotted; cT 1 =4) 
regimes. In the BEC and unitarity regime, we observe an 
almost linear increase of A^ with A - neglecting effects of 
a small anomalous dimension, the universal scaling rela- 
tions indeed suggest A^/A — > const. The linear increase 
of \ v with A extends within a substantial range in tem- 
perature. On the BCS side, in contrast, the universal 
critical region strongly shrinks. 

In Fig. [TUl we show the slope of A v /A on a logarith- 
mic scale. For the BEC and unitarity regime, we can 
numerically read off the slope, A^ ~ A 1 *, £ = 0.98. For 
the unitarity regime (c _1 = 0, dashed dotted line), the 
universal behavior seems to be reached very quickly, typ- 
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FIG. 9: Boson self-interaction A ¥ ,/A ¥ ,(T = 0) as a function 
of the gap parameter A(T)/A(T = 0). The corresponding 
temperature range is < T < T c , where T c corresponds to 
the origin of the plot. Solid line: BCS regime, c _1 = —2. 
Dot-dashed line: crossover regime, c _1 = 0. Dotted line: 
BEC regime, c" 1 = 4. 



d]nA<p(T)/dlnA(T) 
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ln[A(T)/A(0)] 

FIG. 10: Logarithmic slope of \ v , given by dlnA^/dlnA, 
as a function of the gap ln[A(T)/A(T = 0)]. Here T c is 
approached for large negative values of In A. Solid line: BCS 
regime, c _1 = —2. Dot-dashed line: crossover regime, c _1 = 
0. Dotted line: BEC regime, c -1 = 4. 



ically already for A(T) sa A(0)/2. This is related to 
the presence of a strong coupling which induces a fast 
approach to the infrared fixed point. The universal re- 
gion is smaller on the BEC side (c _1 = 4, dotted line). In 
contrast, the approach to the universal fixed point is very 
slow in the BCS regime (c _1 = —2, solid line). In con- 
sequence, the temperature interval around T c for which 
universal behavior occurs is found to be very narrow. 



VII. ATOM DENSITY 

Since we use the grand canonical formalism, the parti- 
cle number n is fixed indirectly by the chemical potential 



p which is a parameter of our microscopic model in Eq. 
([1]) similar as the temperature T or the detuning v. One 
of the dimensionful parameters can be used to set the 
scale of the problem, such that actual computations only 
determine dimensionless combinations of observables as 
functions of dimensionless combinations of parameters. 
To compare our results to experiment or other methods 
one might in principle consider dimensionless quantities 
that involve p, for example T c /p. For example, with the 
truncation described in Sect. [H] we find at unitarity the 
ratio T c / fi = 0.44 while the improved truncation in Ref. 
[23| yields T c //i = 0.39. This may be compared with the 
results from Monte-Carlo calculations T c /p = 0.31 H, 
T c /p = 0.35 @, from 1 /TV-expansion T c /p = 0.232 pj] 
or from 2-PI methods T c /p = 0.41 pf . 

However, the rescaling with respect to p is not op- 
timal for a comparison to experiments. In contrast to 
the particle density n, the chemical potential p is not 
directly accessible experimentally. Dimensionless ratios 
that directly involve the particle number such as T c /Tp 
with T F = (3?r 2 n) 2 / 3 can be measured more directly than 
T c /p. 

For strongly interacting particles the dependence n(p) 
is non-trivial to obtain. This is in contrast to weak inter- 
actions where n(p) can be estimated by the correspond- 
ing formula for a non-interacting gas. Part of the un- 
certainity in the prediction of a dimensionless ratio that 
involves the density arises therefore from the determina- 
tion of the function n(p). In this paper, we use a flow 
equation for a scale-dependent generalization of the den- 
sity, 



d^k = -d k -^-U k , 

Op 



(73) 



in order to calculate the density n = nk\k=o- The fi de- 
pendence of Uk is approximated by an expansion in p and 
li as shown in Eq. (|2"0"]) . For the determination of n, we 
need only the terms linear in 8p, with fj. = ^q+S/j,. Besides 
the term linear in 5/j, and p — po, one might add another 
term which is quadratic in p — pq. We have checked nu- 
merically that the inclusion of such a term has only a 
very small influence on the quantitative determination of 
n and neglect it for our numerical results. 

In this work, we neglect a possible renormalization of 
the fermion propagator G.0. This implies also that the 
dependence on p is the same for all cutoff scales k. In 
particular, the effect of fluctuations on the size of the 
Fermi sphere is not taken into account such that it is 
of radius p 1 / 2 an all scales. Beyond our approximation, 
we expect that the effect of a fc-dependent Fermi sphere 
changes somewhat the dependence of the density n on 
the chemical potential p. 

To illustrate the importance of the density determina- 
tion, we compare our result for T c /T-p with a calculation 
where a much simpler (and naive) method is used to es- 
timate the density. It is motivated physically by assum- 
ing that we have simply coexisting pointlikc atoms and 
dimers, which arc coupled via a common chemical po- 
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tcntial [16| - in this picture, strong correlation effects are 
completely omitted for the density determination, and we 
can get an impression of the sensitivity of our results with 
respect to such correlation effects. At the critical tem- 
perature the superfluid and condensate densities vanish 
and no occupation number density for a single momen- 
tum mode is expected to be of the order n. In addition, 
the energy gap for the bosons just vanishes, mt = 0. Ne- 
glecting the effect of interactions, one could use a naive 
estimate for the density at T = T c , 



d 3 p 
(2tt) 3 



Tc/Tf 



-|»)/T C + 1 e pV(2T c ) _ l 



(74) 

The first term gives the contribution of the fermions 
where the factor 2 counts the spin-states. As stated 
above, the second contribution comes from the bosons 
and the factor 2 reflects now that every boson consists of 
two fermions. Our flow equation method yields a result 
for n that is roughly a factor 1/2 smaller than n na ive- 
Correspondingly, the ratio T C /T F is larger. For a com- 
parison of both methods, we plot T c /Tp as a function 
of cT l in Fig. [TT] At the unitarity point c _1 = 0, the 
formula in Eq. ^ leads to T C /T F = 0.171. 

As a side remark, we note that the approximation 
made to obtain Eq. ([71)1 becomes valid for extremely 
narrow resonances with vanishing Yukawa interaction 
hm — > fl(| . Of course, the strongly interacting problem 
we are interested in here is the opposite limit of broad 
resonances h v — > oo - here, we only try to estimate the 
uncertainty emerging from the neglect of correlation ef- 
fects on the density. 

At the phase transition we have two independent (di- 
mensionless) ratios T c /Tp and n/Tp. For the unitarity 
limit, c _1 = we display these ratios for different ap- 
proximations and compare them to Monte-Carlo results 
in table |TJ The lack of agreement may indicate that the 
truncation should be improved further, and we discuss 
various directions in the conclusions. Another insuffi- 
ciency of the present truncation becomes visible in the 
behaviour of the superfluid density for T — > 0, as shown 
in Fig. 1121 The ratio ns/n should remain smaller than 
one, which is not obeyed by our current truncation on the 
BCS side. At zero temperature Galilean invariance im- 
plies n = ns throughout the whole BCS-BEC crossover. 
We expect that the problem can be solved by including 
the p dependence of the wavefunction renormalization A v 
where the dependence on p was neglected in our trunca- 
tion. 



VIII. CONCLUSION 

In this paper, we have given a rather detailed ac- 
count of the use of the flowing action for ultracold 
fcrmionic atoms. Our method allows for a unified ap- 
proach for all temperatures and densities - including the 
vacuum for T = n = - as well as arbitrary scatter- 
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FIG. 11: (Color online) Dimensionless critical temperature 
T c /Tf as a function of the inverse concentration c _1 = 
(afe) - . We show the result obtained in the truncation de- 
scibed in Sec. [II] with the flow equation for the density (solid 
line). In addition we show the result obtained using the sim- 
pler (naive) estimate for the density in Eq. (|T4[) (long dashed 
line). The dotted line gives the critical temperature of a free 
BEC while the short-dashed line is the expected correction 
for an interacting Bose gas in Eq. (17211. with the Monte-Carlo 
result k — 1.3 and with auja = 0.6 [441 ]. 





T c /Tf 


fi/T F 


Tc/m 


Truncation Sect. [IT] 


0.276 


0.63 


0.44 


Truncation including 


0.245 


0.63 


0.39 


particle-hole flue. [23] 








Density according 


0.171 


0.39 


0.44 


to Eq. (J74J) 








Monte-Carlo [7] 


0.152(7) 0.493(17) 


0.31 


Monte-Carlo [6] 


0.15(1) 


0.43(1) 


0.35 



TABLE I: Critical temperature T c , critical density n c (ex- 
pressed through the Fermi temperature Tp = (37r 2 n c ) 2 '' 3 ), 
and chemical potential /j, in the "unitarity limit" of infinite 
scattering length c _1 = 0. 



ing length. The overall picture of the phase transition 
to the low-temperature superfluid phase and the BCS- 
BEC crossover as a function of the scattering length or 
concentration c = akv is rather satisfactory. In partic- 
ular, the dependence of the critical temperature T c on 
the scattering length yields the effects of interactions be- 
tween composite molecules on the BEC side, as well as 
the correct behaviour on the BCS side including the ef- 
fect of particle-hole fluctuations. The largest uncertain- 
ties occur in the "unitarity limit" c -1 = of very large 
scattering length. In this region in parameter space, the 
fluctuations are strongest and possible deviations of the 
effective action from our rather simple truncation are ex- 
pected to be largest. As far as the critical temperature 
for c _1 = is concerned, the comparison with the Monte- 
Carlo results in table [J may give an idea of the uncertain- 
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FIG. 12: Superfluid fraction ns/n = 2po/n as a function of 
temperature T/Tp. We show the curves obtained on the BCS 
side with c _1 = —1 (dashed line), in the unitarity regime 
with c _1 = (solid line) and on the BEC side with c _1 = 1 
(dashed-dotted line) and c _1 = 2 (dotted line). While the 
superfluid fraction approaches 1 for small temperatures on 
the BEC side, this is not the case on the BSC side and in the 
unitarity regime. This feature is an artefact of our present 
truncation. 



pendence of U on T and /x, as specific heat, compress- 
ibility, or sound velocities. This has been demonstrated 
for bosonic atoms in [29[ and the same methods could 
be taken over to the BEC-BCS crossover system. Fur- 
ther generalizations concern systems with three or more 



ent fermion species 
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species of fermions 40 1 or unequal abundances of diffcr- 



|61| . We hope that the generality 



and flexibility of our treatment will lead to a unified treat- 
ment of many interesting features encountered in these 
systems, and simultaneously permit quantitative accu- 
racy of the theoretical estimates, as needed for precision 
measurements. 
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Appendix A: Units 



ties. 

The phase diagram can be computed from the depen- 
dence of the effective boson potential U(tp) on the tem- 
perature and density or T and fi. The exact flow equa- 
tion (J7|) for the flowing potential Uk(tp) involves the ex- 
act propagators of fermions and composite bosons in the 
presence of a constant background field tp - the conden- 
sate corresponds to the value of tpo for which the mini- 
mum of U{tp) occurs. The present truncation (Eqs. ©, 
(|15l) ) of the propagators remains rather crude and we ex- 
pect a considerable quantitative improvement for trunca- 
tions with a more detailed resolution of the dependence 
of the propagators on momentum and frequency, as well 
as their dependence on tp. On the other hand, increased 
precision also needs an increased accuracy in the parame- 
ters used in a given truncation of the propagators. They 
are determined by exact flow equations for the inverse 
propagators which involve, in turn, the cubic and quartic 
vertices. In this respect the accuracy may be enhanced 
by a more complicated truncation of the momentum and 
frequency dependence of the effective four-fermion inter- 
action, beyond the molecule-exchange channel on which 
we have concentrated in this paper. In this respect, meth- 
ods of rebosonization which closely keep track of the re- 
lation between the composite field tp and its fermionic 
constituents [Bll ] could lead to an impoved treatment of 
particle-hole fluctuations for large scattering length. 

The effort of improving the truncation leads not only to 
an increased accuracy of the phase diagram. Simultane- 
ously, many detailed features of the propagators become 
available. We have demonstrated this briefly by our com- 
putation of the correlation length. Furthermore, many 
thermodynamic quantities can be computed from the de- 



In this work, our units are chosen such that h = A;r = 
2AI = 1, where M is the mass of the fermionic atoms. 
Setting h = 1 implies that space coordinates and inverse 
momenta are measured in terms of the same dimensionful 
scale. The choice /cb = 1 similarly establishes an equiv- 
alent measurement of energy and temperature in terms 
of the same scale. Setting 2M = 1 finally relates the 
measurement scales of energy and temperature on the 
one hand to squared momenta or inverse squared coor- 
dinates on the other hand. (The meaning of the latter 
choice is similar to setting the velocity of light c = 1 in 
relativistic theories). Only one unit remains free, say the 
length I, with dimensions [t] = I 2 , [q\ = [T] = /~ 2 , 
[/i] = r 2 , and [n] = l~ 3 . 

In practice, these units are obtained by rescaling the 
coordinates x/h—> x for H = 1, the temperature k^T — > 
T for fee = 1, and the time t/2M -> r for 2 M = 1. 
The last step goes along with rescaling also all interac- 
tion terms of the action: 2MS- m t — > Si n t- All dimen- 
sionful quantities are then measured in terms of only one 
dimensionful scale which we set in terms of the density 
n of the system or, equivalently, in terms of the Fermi 
momentum 

fc F = (S^n) 1 ? 3 . (Al) 

Typical values of the Fermi momentum occurring in 
trapped Fermi-gas experiments are of the order kp ~ 
0(1 . . . lOeV/c) in units of electron Volt/c (The light ve- 
locity c occurs here to convert the energy unit eV into a 
momentum). With our choice of units, = leV/c cor- 
responds to a density n ~ 4.4 x 10 12 cm -3 for which we 
have used that leV/c ~ 5.07 x 10 6 m _1 in our units with 
H = 1. A typical length scale such as the Bohr radius 
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as = 5.29177 x 10 11 m then reads in inverse momentum 
units a B = 2.6817 x KT^eV/c)- 1 . 

Energy- or temperature-like quantities can most con- 
veniently be measured in terms of the Fermi energy or 
Fermi temperature which are equivalent in our units, 



As 1 = 2M 



k 2 

— rr "'F 7 2 

2M F 



M 



(A2) 



J cv/c 2 e ^ / c ' nxm S the atom mass M in 
units ofeV/c 2 (e.g., Me Li ~ 5.6xl0 9 cV/c 2 for Lithium-6) 
corresponds to expressing the velocity of light c in units 
of VcV. This establishes cV units for our cncrgy-likc 
quantities such as the Fermi energy. 

As shown in the main text, the scattering length a 
remains as the only physical parameter of the interacting 
Fermi gas in the broad-resonance limit. Together with 
the temperature, the relevant dimcnsionless quantities 
characterizing the phase diagram are 



T 
TV 



T 
fcp 



(A3) 



where c (not to be confused with the velocity of light 
here) can be viewed as a concentration parameter, as 
it measures the inverse scattering length in units of the 
interparticlc spacing ~ kp 1 . 

We use our result for the critical temperature at uni- 
tarity, T c /Tp ~ 0.276 at akp = 0, to illustrate how to 
convert our results into standard laboratory units, us- 
ing Lithium-6 as an example. First, we assume an atom 
density of n = 10 13 cm -3 which implies kp ~ 1.3cV/c. 
The resulting critical temperature reads T c ~ 0.276Tp ~ 
0.477(eV/c) 2 . In our units, where the Lithium-6 mass 
obeys 2M = 1, the velocity of light can be expressed as 

c = ^2-^j^cV ~ 1.06 x 10 5 VcV. Inserting this into 

the critical temperature yields T c ~ 4.26 x 10~ 11 eV ~ 
0.49^K. 



Appendix B: Symmetries 

Symmetries of the microscopic action S are also sym- 
metries of the effective action T, provided the functional 
measure does not give rise to an anomaly and that the in- 
frared cutoff respects the symmetry. The symmetries of 
r severely restrict its possible form. They are therefore 
helpful in devising useful truncations. All truncations 
should be consistent with the symmetries. 

In our case we have first the symmetries of translations 
in space and time and the rotation symmetry in space. 
Further, a continuous U(l) symetry reflects the conserved 
number of atoms. 



02 ' 



e">2, 



02 



V 2 * 



This symmetry is conserved by the cutoff and easily re- 
alized in truncations by including in V only terms with 
total atom number zero. (Here ip and ip have atom num- 
ber one and two, and ip* , ip* the corresponding negative 
numbers.) The action {1} also shows an SU(2)-symmetry 
of "spin rotations" between the two fcrmion compo- 
nents. This is again easily implemented. Further discrete 
symmetries are parity reflections (i/j(t,x) —¥ 0(t, — x), 
ip(r,x) — > ip{T, —x)) as well as cuclidcan time reflection 

i/>i(t,x)-> 1 (-r,f), ip{{r,x) -> — 0i(-r,f), 
ip 2 {T,x)^ -tp*(-T,x), ip£(t,a!)->ip2{-T,£), 
<p(t,x)-> <p*(-T,x), <P*(t,x) -> <p(-r,x). (B2) 

We note that ip a and 0* are separate Grassmann vari- 
ables for which the structure of complex conjugation is 
not defined a priori. The transformation (|B2[) is actually 
not compatible with the usual complex structure. (The 
action (Q} is also invariant under ip(r,x) tp*(—T,x), 
ip a (r, x) o ip a (— r, x) if all Grassmann variables are com- 
pletely reordered.) Time reflection symmetry does not 
forbid terms with an even number of r-derivatives. 

For zero temperature no rest frame is singled out 
a priori and we expect invariance under Galilei- 
transformations. The implications of Galilean invariance 
can be understood best if we perform an analytic con- 
tinuation from euclidean time to real time r = it. The 
microscopic action |T]) becomes then 

S[tp, ip)=i j dt j d 3 x(ip*< {-id t - A - n) 



(B3) 



-tp* ( —id t - 2 A + v - 2(i) (p 



-h v {tp* i\) T tip - tptp'eip*) 



Galilean boost transformations (with a boost velocity 2q) 
act on the fields as 

ip(t,x)^ip'(t,x) = e~ l ^ 2t -^iP(t,x-2qt) 



ip(t, x) — > tp'(t, x) = e 



-2i(q 2 t-qx) 



tp(t,x-2qt). (B4) 



While the invariance of the Yukawa term under this 
transformation is obvious, its realization for the kinetic 
term is more involved. Performing the transformation 
explicitly one finds 

ip^Aip -> ifi Aip - cfip^ip + 2iqip^Vip 

tp^ Atp — > ip^ Aip — Aq^ip^ tp + Aiqtp^Vtp 

tp'idtip — > ip^idtip + (fip^ip ~ 2iqip jf Vip 

ip'id t tp —> tp'idttp + 2(fip' { 'ip — 2iqip !f Vip 1 (B5) 

such that indeed the combinations 



for the fcrmions and 



ip* -> e- 2la ip* 



(Bl) 



id t 



id t 



(B6) 
(B7) 
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for the bosons lead to this invariance. The conserved 
Noether charge of Galilean symmetry is the center of 
mass momentum. 

In the case of a bosonic condensate tp(x, t) = tpo the 
Galilean symmetry is spontaneously broken. Now the 
condensate determines the rest frame. The spontaneous 
symmetry breaking of Galilean symmetry comes in pair 
with spontaneous symmetry breaking of the U(l) symme- 
try associated to atom number. It can be associated with 
superfluidity. We stress that for T > Galilean symme- 
try is explicitly broken since the heat bath singles out 
a rest frame. For low temperatures truncations consis- 
tent with Galilean transformations (and therefore always 
containing the combinations (|B6|) and (|B7|) remain use- 
ful, however, since symmetry breaking terms vanish for 
T — > 0. 

Finally, we observe another useful invariance of the ac- 
tion for T = 0. It arises if we consider p as a source 
term for the fermion and boson bilinears and consider 
transformations that also change p. (These type of in- 
variances are not symmetries of the action in the stan- 
dard sense, since "parameters" of the action are changed 
by the transformation.) Let us extend Eq. (|B3|) to a 
t-dependent source p{t). There is now a semilocal (7(1) 
invariance of the form 

, 2ia(t) * , —2ia(t) * 

cp — > e <p, <p — > e y '(p 

p -> p + d t a. (B8) 

This holds, since the combinations {—idt — p) and {—idt — 
2p) act as covariant derivatives for the fermions and 
bosons, respectively. 

Combining semi-local U{1) symmetry and Galilean 
symmetry at T = we find, that the derivative operators 
and the chemical potential are combined to the operator 

D4, = {-idt - A - a*) (B9) 
for the fermions and 

D v = {-id t - — -2fi) (BIO) 

for the bosons. In addition to powers of this operator, 
only spatial derivative of terms that are invariant under 
U{1) transformations like e.g. pAp with p = ip*<p may 
appear. The above symmetry transformations act lin- 
early on the fields and so also the effective action V is 
invariant. Also there, only powers of the operators 
and D v may act on the fields ip and ip. 

In the presence of spontaneous symmetry breaking by 
a nonvanishing order parameter p = ip*ip > the possi- 
bilities for combining space and time derivatives get more 
complex. Effectively, terms quadratic in dt need not be 
matched by terms quadratic in A or /i. For a more de- 
tailed discussion see [28| . Let us finally mention, that 
there is an additional conformal invariance at the uni- 
tary point in the crossover where the fermion scattering 



length diverges c = akp = oo [57j • The only length scale 
is then set by the chemical potential p or by the particle 
number density n. This gives further constraints to the 
effective action at that point in the phase diagram. 

Appendix C: Cutoff functions 

The computations in the present work are done with 
regulators which are optimised w.r.t. the derivative ex- 
pansion, see e.g. p58l - [60| . Such a choice stabilises the 
truncation and minimizes the systematical error for the 
computed obscrvablcs. For the scalar field we choose the 
Litim regulator [58| 

i?^V) = A v k 2 r kv {y), V=^, 

r kv {y) = {l-y)6{l-y). (CI) 

Combining the g^-dependent term in the inverse prop- 
agator G^ 1 with the cutoff term one finds for P v = 

q s /2 + R[ v] " jA v the simple behavior P v = cp /2 for 
q 2 > 2k 2 , P v = k 2 for q 2 < 2k 2 . In other words, our 
regulator leaves the propagation of modes with spatial 
momenta bigger than the cut-off scale, q 2 > 2k 2 , un- 
changed. For modes with spatial momenta smaller than 
the cut-off scale, q 2 < 2k 2 , it effectively introduces an 
infrared mass term k 2 . 

For the fermionic field we choose an analogue of the 
scalar regulator (|C1|) . Here, however, the propagation of 
modes close to the Fermi surface is suppressed, 

i*«V) = k 2 r^{z), 
r ki> {z) = {sign{z)-z)e{l-\z\), (C2) 

where 

z = {q 2 -p)/k 2 . (C3) 

This regulator leaves the propagation of modes with spa- 
tial momenta farther away from the Fermi surface than 
the cut-off scale, q 2 — p > k 2 , unchanged but leads to a 
flat (inverse) propagator k 2 for modes with spatial mo- 
menta closer to the Fermi surface than the cut-off scale 
q 2 — p < k 2 . This effectively introduces an infrared mass 
term proportional to k 2 . In summary this takes into ac- 
count the fact that fermionic fluctuations have to be cut- 
off relative to the Fermi surface, and hence introduces 
to leading order the same cut-off scales for the fermionic 
and bosonic fluctuations. The fermionic regulator in the 
presence of a Fermi surface is displayed in Fig. [T3l 

Appendix D: Flow of the effective potential 

In this App. we derive the flow equation for the ef- 
fective potential U{p,p). We start from the exact flow 
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► 9(1 



t 6(1 - |z|) = 



S-8(1 - |z|) = 1 



sign(z) = - 1 



sign(z) = +1 



FIG. 13: The fermionic regulator in the presence ol a Fermi 
surface p > 0. 



equation for the average action (or flowing action) [27| 
d k T k = ^STr(lf > + R^duRu. (Dl) 



This is evaluated for a bosonic field (p that is constant 
in space and in (Matsubara-) time, tp(x, r) = tp. (The 
expectation value of the fermionic field ip vanishes - for 
vanishing fermion source - due to its Grassmann prop- 
erty.) Using the truncation in Eq. ([13)1 we find the flow 
equation for the potential 



d k U k \ p = \ 



[(G kv )n + {G kip ) 2 2\ A^dk {A v k 2 r kv ) 

{Gkijj) 31 — {G k ^,)i2\ 



[(Gktp)l3 + (G k ^)2i 

x5 fc (k 2 r k ip) \, 



(D2) 



which is evaluated for fixed p = <p*p. The dimensionlcss 
function r klp depends on y = q 2 /(2k 2 ) while rkiji depends 
on z = (q 2 — po)/k 2 . The regularized propagators G ktp 
and G k ^p that appear in Eq. |D2]) arc modified from Eqs. 
(fTTT) (|15|) by the presence of the ultraviolet regulator 



G k l 



= G~ 




G kt f, — G y 



k 2 r k ^ 





k 2 r k(p J 

~k 2 r k4 , 




(D3) 



We can now perform the summation over the Matsubara 
frequencies go — 2ixTn for the bosons and go = 2irT(n + 
1 /2) for the fermions. The integration over <f is performed 
quite generally in d spatial dimensions. The result can 
be expressed in terms of the dimensionless variables 



Wi = 

Mo 
k 2 ' 



Ek 
k 2 ' 

Ap 



w 2 



2pU" 



V ~ Po 
k 2 ' 



k 2 
f 



T 
k 2 ' 1 



w 3 



k± ' 
= ^ (D4) 

sl<n 



and the anomalous dimension 

kd k A v 



A. 



(D5) 



A term involving the anomalous dimension appears if we 
evaluate the flow for fixed p = A v p instead of fixed p. 
The flow of the effective potential reads now 



k& k U k = r) A oU k + 8V2 



k d+2 



Vd, 



,k d+2 v d , 
8 — l{p)s 



dSif 
(o) 

F ' 



1 - 



d + 2 



.(0) 



(D6) 



Here, the coefficient Vd is proportional to the surface of 
the (i-dimensional unit sphere, which is (2tt) 4vd, with 



1/(8tt 2 ). 



2 d + 1 ir d / 2 T(d/2). In particular one has v 3 



The threshold functions s 



(0) 



,(0) 



(w 1 ,w 2 ,T,S< p ,r]A v ) 
and Sp ^ = s F °^ (w 3 , fx, A/i, T) as well as the function l(p) 
used in Eq. (|D6[) depend on the choice of the infrared 
regulator functions r kip and r k ip- They describe the de- 
coupling of modes when the effective "masses" Wj or —jj, 
get large. The threshold functions for the bosonic fluc- 
tuations reads 



1 - 



d + 2 



AO) 



= dldyyi 
o 

i(u>l + W 2 ) 



{r kv - yr kip - T] Aip r kv ) 



' kip 



y/wi + y + r kip ^/w 2 + y + r kip 



- + N B {yJwi + y + r rvx /w 2 + y + r k(p /S v ) 



and involves the Bose function 



N B (e) = 



1 



e^ T - 1 

and the optimized cutoff function 



r kv (y) = (l-y)e(l-y). 



(D7) 



(D8) 



(D9) 



We normalize the threshold function such that Sg = 1 
for w\ = w 2 = T = 0. For our cutoff one finds the 
particular simple expression 



,(0) 



1 



Wl 



1 + w 2 



1 



Wl 



- + N B Wl + w iy /l + W2/S v ) 



(D10) 



The threshold function for the fermionic fluctuations 
is obtained similar. For a generic cutoff that addresses 
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the spatial momentum, it reads 

/oo 
dz (z + p) (r fe ^ - zrj^) 

(z + r fe ,0 - A/2) 



+ (z + r fc ,/, - A/i) 
^ - A^F w 3 + {z + r hl p - Ap,) 2 ^ 



Here we employ the Fermi function 

1 



N F (e) 



(Dll) 



(D12) 



Note that for a generic cutoff the right hand side of equa- 
tion (|D11[) does not necessarily factorize. In that case 

one might work with a threshold function 4°^ that also 
depends on fl. For our cutoff 



r kM> = (sign(z) - z)6(l - \z\). 



(D13) 



one has Sp = 1 for w 3 = p, = A/i 



and therefore A/i 
the simple form 

- 



0. For /i = /io 
the threshold function assumes 



2 - N F {VT 



w 3 



while 



m = 9& + m + if 2 en* i)(a - i) d/2 . (Du) 

In the limit T = T/k 2 — >• the thermal contributions to 
the flow of Uk vanish, N-q = N F = 0. 

Taking a derivative with respect to p on both sides of 
Eq. (|D6[) we obtain 



kd k U' k = VA (U' k + pU^)-8V2 



k d v d 



dS, 



1 - 



d+2 



U£s%> (S > + $U'{ + 2pU£>)s 



k °B 
k d-2 



(3)^(0,1) 



B 



V d,2 U ~\ C 1 ) 



(D15) 



Here we introduce the derivatives of the threshold func- 
tions 

d 



(0,1) 



.,(D 



dwi 
d 

d 
dw 3 



(0) 
'B 



,(0) 



(0) 



(D16) 



We may divide these into contributions from quantum 
and thermal fluctuations 



,(1.0) 


= (W2 


i (i,o) 


+ s (1 ' 0) 




,(o,i) 

& B 


= (w 2 


\ (0,1) 
~ w l) s B,Q 


' b B,T > 




F 


~~ *F,Q 


+ s (1) 




(D17) 



For T — > the thermal contribution vanishes 4t' = 

s b't = s ft = 0- We extracted a factor (tu 2 — u>i) from 

the threshold functions 4q and Sgg to make explicit 
that these contributions vanish for wi = 102, which holds 
for /) = 0. 

For our choice of the regulator functions r k ,iji and r k tp 
we find the explicit expressions 



,(l,o) 



B,Q 



,(0,1) 

S (1) 
F.Q 



4(l+u>i) 3 / 2 (l + u;2) 1/2 ' 
1 

_ 4(l + u;i) 1 /2(l + U ; 2 )3/2' 
1 



2(1 + w 3 ) 3 / 2 ' 



For T > the thermal fluctuations lead to the additional 
contributions from the bosons 

4,t = 2( W2 - Wl ) 4,q ) N B (v/(l + Wl )(l + w 2 )/S v ) 

+ S -Q;ra;^( v(1+Wl)(1+ ^)' 

4°T = 2( W2 - t»i) 4 1 ,Q ( V(l + ^l)(l + ^2)/S v ) 



+ 4 0) Q 7=| 7V B + *>i)(l + W2 )/S v ) . 

(D18) 

Here, we use the derivative of the Bosc function 



N B (e) = ^iV(e). 



(D19) 



Similarly, the fermionic part of the thermal contribution 
reads 



S (D 

F,T 



„(o) 1 



with the derivative of the Fermi function 

N F (e) = §-N F (e). 
ae 



K (VT+Ws) , (D20) 



(D21) 



Appendix E: Gradient Coefficient and Wave 
Function Renromalization 

We compute the flow equations for the kinetic coef- 
ficients A v , Z v for the bosons according to the follow- 
ing procedure: First, we obtain the frequency and mo- 
mentum dependent one-loop correction to the inverse bo- 
son propagator in the presence of the infrared cutoff R k . 
The flow equations for the boson propagator are then 
obtained by taking the appropriate cutoff derivative. Fi- 
nally, we project onto the flow equations of A v and Z v by 
taking derivatives with respect to the momentum q 2 and 
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the frequency qo, respectively. This procedure is equiv- 
alent with the derivation of the exact flow equation for 
the inverse propagator [27j , [32j , and a truncation where 
only the couplings appearing in Eq. (|13[) are taken into 
account. 

We work in momentum space with a basis of bare real 
field fluctuations 8tp\,8(p2, which relate to the complex 
bosonic fields in momentum space as 

8<p{Q) = (8(p 1 (Q) + iS(p2(-Q))/V2, (Ef) 
8(p*{Q) = (8<pi(-Q) - iS<p 2 (Q))/V2. 

The one-loop correction to the inverse boson propagator 
AP is then defined with 



(AP v ) ab (q)8(q 



S 2 AT 



8<p a (q)8(p b (q') 



(E2) 



Here Ar is the one-loop correction to the effective action. 
It is evaluated for a constant background field (p\ = \J2p§ 
if the minimum of Uk occurs for po(k) 0, while (p\ = 
in the symmetric regime where po{k) = 0. Then the wave 
function renormalization and gradient coefficient at one 
loop order are obtained as 



AZ ir . 



■(AP v ) 12 (g ,0) 



90=0' 



(E3) 



A ^ = 2 ^(Ai>,) 22 (o,g-)| f=0 . 



Within our truncation, (AP V ) 12 and (AP ip ) 2 2 receive 
contributions from fermionic and from bosonic fluctua- 
tions, 



(AP„)ia(gt,,0) = (Ap( F )) 12 (g ,0) + (Ap( B )) 12 ( (Z o,0), (E4) 

(A^ F) )u(?o,0) = q hl [ { ( f ^ 

v J p det F (p , p) det F (p + qo,P) 

/a5(B)n / n \ 7 \2 - f 3a 22 (p) - an(p) 
{APS, J )i 2 (<7o,0) = Z v q \ p 1 



p det B (po j P) dct B (po + 3o , p) ' 
(AP ¥ ,) 22 (0,g) = (APf)) 22 (0,<?) + (Ap( B )) 22 (0,g), 

v Jp det F (po , p) det F (po , P + q) 



, A p(B), m -* T2- f 2(Z y p ) 2 + an(p + g)a 22 (p) + aii(p) a 22 (p + g) 

( A Pp J22(U,(?J - -A po / j— 7 it-; i ; _ , -> 

J P det B (po,Pjdet B (po,P + «} 



Here we use the abbreviations 

f(p) = p 2 -p + R«>\p 2 ), (E5) 

det F (u,p) = p 2 + f 2 (p) + h*po, 

a n (p) = A v p 2 /2 + 2~X v p + R^\f), 

a 2 2(p) = A v p 2 /2 + +R<j?\f), 

dets(w,p) = (Z v p ) 2 + an(p)a 22 (p"). 



(We work with the boson propagator parameterized ap- 
propriately for the regime with U(l) symmetry breaking. 
In the symmetric regime the contributions from bosonic 
fluctuations vanish.) To obtain the flow equation for the 
boson propagator (P v )i 2 , we take the cutoff derivative of 
the expressions in Eq. (1E5|) 
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(d fc iy i2 (q ,0) = (d k P^) 12 (q ,0) + (d k P^y) 12 (q ,0), 

qo - hl f f fctfV) 



dct F (po , p) det F (p + qo,p) 
2/(p) 2 d k R ( f (?) 2f(p) 2 d k R^ {f) 



det F (p ,p) det F (p + q ,p) det F (po , p") det F (p + qo,p) 
(3 fc P< B) )i 2 (<Zo,0) = 4(AP( B >) 12 (<z ,0) 



det B (po,p)dct B (p + qo,P) 



[3ag2(g) ~ aii(p)][an(p) + a 2 2(p)]^4 y) (P g ) 

det B (p , P) det B (po + <7o , p) 
[3a 22 (p) -aii(p)][an(p) +a22(pl]^fij y) (p 2 ) 
det B (p , p) det B (p + <?o , p) 



(E6) 



Here the symbol 9& denotes a formal derivative that hits only the cutoff term R k . Similarly, we obtain the flow 
equation for (P v ) 22 as a function of the spatial momentum q 

(d k P v )22(0,q) = (d k P^) 22 (0,q) + (d k P^) 22 (0,q), 
(d k PP) 22 (0,q) = d k (APP) 22 (0,p) 

= _- h 2 j[ f{p + <t)d k R { ^\?) 

v Jp \ dct F (p , p) det F (p , p + q) 

2b§ + /(p)/(p+ g) + hiMmdkR^m 

det F (p , p) det F (p , p + (?) 



detF (po,p — q) dct F (p , p) 

det F (p ,p- q)dct 2 F (p ,p) 



(d k P^) 22 (0,q) = d k (ApW) 22 (0,q) 



■ [ail(p + g) + a 22 (p + q}]d k R^i?) 



det s (p , p) det B (po , P + q) 

[2(Z y p ) 2 +Qii(p + g)a22(p) + «ii(p)a22(p + g)][an(p) + a22(p)]d k R i k ) jp 2 ) 
det 2 B (p , p) det B (po , P + q) 

[ ai l(p - 3) + a 22 (p - 



(E7) 



dets (po,p — q) dets (po , p) 
[2(Z y p ) 2 + aii(p)a 22 (p - g) + an(p - g>22(p)][an(p) + a 22 (p)]d k R k v) (p 2 ) 
det B (po,p-q) det 2 B (p , p) 

I 

In Eq. (|E7[) we performed a shift in the integration vari- Z v and ^ by taking approriatc derivatives 
able p — s- p — <f at several places. From Eqs. (|E6j) . (|E7j) g 

one obtains the flow equation for the kinetic coeffcients d k Z v = — — (d k P ip )i 2 (qo,0)\ q _ Q , 

d k A v = -^(d k P v ) 22 (0,q=(q,0,0))\ ?=o - (E8) 
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After this projection and using the cutoff functions in 
Eqs. (|C1[) . (|C2[) we can perfom the Matsubara summa- 
tion over the frequency po and the integration over the 
momentum p. We also switch from the bare couplings to 
the renormalized ones according to the discussion below 
Eq. (fT3|) ; the latter are obtained from rescaling with the 



flowing gradient coefficient A v . The effect of the gradient 
coefficient is then encoded in the anomalous dimension 
Eq. (|D5j) . We display the flow equations in the symme- 
try broken phase, and separate into the fermionic and 
the bosonic contribution, 



(F) . (B) 



„(F) 



>)A 



67T 2 fc-y/l + W3 



3/2 



37T 2 fc(f + W 2 ) 3 ^ 2 S 2 



(fi + lf' 2 Q{~p + 1) + (/} - lf' 2 0{p - 1 )] (1 - 2iV F (v / m^) + 2 V / TT^<(VTT^ 3 ")) , 

(E9) 



([1 + 2N B (VT+^/S v ,)]S ip - 2VTT^Ni 3 (VTT^/s v )). 



The Bose and Fermi functions A^b, Np are defined in Eqs. (|D8[) . (|Df 2[) . and the derivative refers to the full argument. 
The other abbreviations are defined in Eq. (|D4[) . The flow equation for the rescaled wave function renormalization 
then reads 



d t S v = p s + r, Av S v , p s = 4 F) + Pf ] , 

= M( fl W3)5/2 ( 2 ^(M) - (A - D^A -!)-(£+ D 3/2 ^(M + D) 

x (27V F (Vl + w 3 )(w3 - 2) + 2(1 + 2V1 + W3^f(v / 1 + w 3 )) 

+2w 2 ^'(\/rT^3") - ^3(1 + 2VTTmN F {VTT^) - 2N$(VT+wZ))), 

4 B) = 3^(^^5/252 ~ 2^Vb(v / T+^2'/^))^ + 2 v / TT^^jV B (v / T+^/^)](2 -W2- 3/4^) 

+w 2 (2 + 5/2w 2 + wl/2)N£(y/T+m/S v )) . (EfO) 



As mentioned above, the bosonic contributions to the 
above flow equations are ~ po , thus vanishing in the sym- 
metric phase. The fermionic contributions in the latter 
case are obtained via the replacement po — > 0, imply- 
ing W3 —> 0. In the zero temperature limit T/k 2 
the thermal contributions to the flow equations vanish, 
formally obtained via the replacement Ab = Af = 0. 

Appendix F: Flow of the Yukawa coupling 

In our truncation Eq. (fT4| the flow of the Yukawa 
coupling is quite simple. In the symmetric regime with 
po = there is no loop contribution such that the flow is 
given by the anomalous dimension only, 

kd k h v = 0, kd k h v = ^r] Av h v . (FI) 

For a nonvanishing order parameter po > there is a 
loop contribution ~ h^X^po from a diagram involving 
both fermions and bosons. This contribution is expected 



to be subleading and we checked numerically that this is 
indeed the case. We have omitted this contribution for 
our numerical result. 

The situation is different for an extended truncation 
where also the effect of particle hole fluctuations is taken 
into account. This was described briefly in Sect. [V] and 
in more detail in [23| . The box-diagram Fig. [7] gives a 
contribution to the flow of the four-fermion interaction 
A,/, even though A^, = was realized on the initial scale 
k = A by virtue of a Hubbard-Stratonovich transforma- 
tion. On the microscopic scale the only contribution to 
the interaction between fermions is given by the exchange 
of a bound state <p, i. e. a term as in Eq. ([2]). Using 
the formalism of rebosonization [48( one can maintain 
this picture also for k < A and adapt the flow equa- 
tion for h v (and possibly m 2 ) such that the effect of the 
box-diagram in Fig. [7] (the particle-hole contribution) is 
taken into account. This setting where the regeneration 
of the four-fermion coupling is suppressed in favor of 
a modified flow equation for the Yukawa coupling h v is 
descibed in more detail in (23j . 
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